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to combustion dynamics, chamber dynamics, engineering model 
formulation, and computer programming, and Mr. K. W. Fertia with 
respect to numerical analysis, computer programming, and checkout. 
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ABSTRACT 


This report summarizes the technical effort conducted during an eleven- 
month program to develop and verify a digital computer model for NASA/ JSC 
which can be used to analyze and predict engine/feed system coupled insta- 
bilities in pressure-fed storable propellant propulsion systems over a 
frequency range of 10 to 1000 Hz. 

The analytical approach tc modeling the feed system hydrodynamics, com- 
bustion dynamics, chamber dynamics, and overall engineering model structure 
is described and the governing equations in each of the technical areas is 
presented. This is followed by a description of the generalized computer 
model in which the specifics of the hydrodynamics, combustion dynamics, and 
chamber dynamics are formulated into discrete subprograms and integrated 
into an overall engineering model structure. 

The operation and capabilities of the engineering model are verified by 
comparing the model's theoretical predictions with experimental data from 
an OMS-type engine with known feed system/engine chugging history. The 
latter data were obtained at White Sands Test Facility on Rocketdyne hard- 
ware during Task XII of the SS/OME Reusable Thrust Chamber Program 
(NAS9- 12802). 

The program is concluded with the successful operation of the engineering 
model on the NASA/ JSC Uni vac 11 10, EXEC-8 computer system and by extensive 
documentation of the model in the form of a computer user's guide and 
final report. 
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SECTION I 


INTRODUCTION AND SUMMARY 

Historically, during the development of pressure-fed propulsion systems, 
feed system/engine coupled instabilities have been frequently encountered. 
Resolution of these problems usually included increasing injector pressure 
drop to decouple the feed system from the combustor, the result being 
substantial system weight penalties. A dynamic computer model would be 
a useful tool in obviating coupled stability problems during the 
ment of the Space Shuttle Orbit Maneuvering System (SS/OMS). A model could 
be used both as a system design tool to optimize component location and 
pressure profile (minimize system weight) and a system development tool co 
define test programs for assessing stability margins of the OMS. 

Models have been constructed to study specific problems on specific engine 
configurations. Under a previous contract, NAS9-10319, a generalized 
propulsion system model for very low frequencies was constructed to provide 
a modular digital computer program for simulating transient operation in 
pressure-fed rocket engine systems (Ref. 1). The use of the program wa* 
demonstrated by modeling the Apollo Ascent and Descent engines and the 
Rocketdyne SE5-5 Propulsion System. 

This document is the final report of an eleven-month program conducted bv 
Rocketdyne to develop and verify an engineering digital computer model for 
the NASA/ JSC which can be used to analyze feed system/engine coupled in- 
stabilities in pressure-fed, storable propellant, propulsion systems over 
a frequency range of 10 to 1000 Hz (frequencies lower than the chamber 
transverse frequencies). The model is sufficiently aener^ so that it may 
be readily applicable to present and future engine and propulsion programs. 
For scaling purposes the baseline conflyuratioi cho;en is the OMS engine. 
The model has been written for use on the NASA/ JSC Ini vac 1110, EXEC-8 
system, and provides NASA a tool which can be used to: 
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1. Conduct preliminary design tradeoff fur feasibility 
studie c prior to propulsion concept selection. 

2. Guide the design of propulMon systems to ensure stability 
at all operating ranges and with minimum penalties. 

3. Guide testing programs by predicting the least stable 
operating regimes thereby reducing e number of stability 
tests required. 

4. Provide stability verification in the event system changes 
are made and hot-fire verification is impractical. 

5. Diagnose problems on existing systems and evaluate potential 
solutions. 

The specific end products of this effort are as follows: 

‘i . A digital generalized computer ,iodel for inves*' tinq 
coupled instability in a Dressure-fed propulsion system. 

The mood is suitable for use in evaluating the stability of 
of the Space Shuttle OMS (Orbital Maneuvering System) in 
its various flight configurations and over its planned 
operating conditio 0 . The structure of the model is modular 
such that it can be adapted to various configurations without 
ma v or modifications. The output format of the model is such 
that the margin of stability is apparent rather than a simp'} 
stable/unstable prediction. Specific frequencies and modes of 
oscillations are also obtainable from the model. 

2. Checkout of the model on the NASA/OSC ( ’ohnson Spacecraft Center) 
Univac 1110, EXEC-8 computer and verification by comparing of a 
WST r uMS enqine and test rig with known feed sys cem/enqlne chug- 
ging or buzzing history. 

3. Complete documentation of the model effort including: 

a. A final report describing the pfforr ir. olved in 
developing the model. 

b. A cc iter manual consistina of 

• ser's section designed for the engineer who desires 
to apply the model to a given engine system. 
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• A programmer's section for elucidation of the programming 
details, i.e., format, subroutining structure, program 
numerical limits, etc. 

• An engineering section explaining the physical implications 
of the mathematics involved, the assumptions, range of 
validity and realistic limitations. 

4. Discussions and recommendations for additional areas of investi- 
gation either for the improvement in the model or in application 
of the model to benefit existing or future propulsion programs. 

One such item is a test plan for the empirical determination of 
appropriate values of the variables describing the injection and 
combustion process. Areas of application of the model include 
analysis of the several OMS flight configurations, investigations 
of any coupled stability phenomena which might be incurred in 
the OMS system to determine their causes and to arrive at the 
most expeditious solutions, and to elucidate the coupled stability 
margin as a function of operating conditions. 

The program was accomplished in three phases: 

I. Model Evaluation anj Formulation 

II. Computer Model Development 

III. Model Verification and Documentation 

A logic flow diagram of the completed effort is shown in Fig. 1 . 

In the initial phase, a detailed assessment of the available techniques 
Tor modeling the propulsion system's hydrodynamics, combustion dvnamics 
and chamber dynamics was conducted. Recommendations for the character- 
ization technique to be used in each of the technical areas were made and 
the dynamic equations were formulated. The results of Phase I were presented 
at a briefing held at NASA/JSC on 15-16 November 1974 and surmarized in 
detail in a Phase I report, "Engineering Model Characterization Evaluation 
Interim Report" (Ref. 2 ). 
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PHASl i 


PHASE II 


PHASE III 



Figure 1. Logic Diagram of OME/Feed System Coupled Stability Program 















In the second phase, the specifics of the hydrodynamics, combustion and 
chamber dynamics were programmed as discrete subprograms and integrated 
into an overall engineering structure to provide a usable, cost efficient 
computer package. Each subprogram was debugged and checkout cases were 
run to ve~ify the computer logic formulation. The overall model was satis- 
factorily executed and reasonable values obtained for the model output 
parameters. A program review meeting was held at NASA/ JSC on 10 April 1975 
to support the results of Phase II. 

In the third phase, the operation and capabilities of the engineering model 
were verified and the model's theoretical predictions were compared with 
experimental data from an OMS-type engine with known feed system/engine chugg- 
ing history. The latter data were obtained at White SandsTest Facility 
(WSTF) on Rocketdyne hardware during Task XII of the Space Shuttle/OME 
Reusable Thrust Chamber Program (NAS9-12802) . The particular sequence of 
data chosen involved a series of integrated thrust chamber tests in which 
the level of chamber pressure was progressively reduced until chugging 
occurred and included the effect of mixture ratio on the minimum pressure 
levels. All pertinent WSTF hot-fire test data to be used durino model 
verification were summarized in an "Engineering Model Verification Plan," 
submitted to NASA/JSC on 7 May 1975. Included in the report were details 
of Rocketdyne's thrust chamber and injector (like doublet No. 1), schematics 
of the fuel and oxidizer configuration, and complete steady-state operating 
data for each of the seven verification analyses which were conducted. 

A total of seven model verification analyses were completed using the Feed 
System Coupled Stability Model based upon the WSTF integrated chamber low 
pressure tests. In each gi-.an test case, the freqc.ncy and stability pre- 
dicted by the engineering model were compared to the experimental data. 

The model was found to be in agreement with the experimental data in pre- 
dicting engine stability or instability in six of the seven analyses. 

In comparing the experimental and calculated frequencies for the unstable 
tests, the calculated unstable frequencies we^e found to be approximately 
16% lower than the experimental values. These discrepancies can be at- 
tributed in part to the assumptions used in modeling the feed system and 
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to the lack of experimentally determined combustion characterization 
parameters which are required as input data to the model. The calcula- 
tions completed thus far have shown that the model is highly sensitive 
to changes in the hydrodynamics system as well as to certain combustion 
parameters (such as the Klystron constant). 

The program was concluded with the conversion of the engineerino model 
from Rocketdyne's IBM 370 computer to the NASA/JSC Univac 1110, EXEC-8 
computer system and successful operation of the engineering model usinq 
the WSTF verification analyses. Model documentation in the form of the 
present final report and a computer user's manual constituted the end 
products of the contract. 

The work performed within all of the foregoing tasks is summarized in 
this document. The presentation of the subject matter is organized as 
a task-by-task description rather than a detailed discussion of the 
computer program. The latter is extensively described in a separate 
companion document, entitled "OME/Feed System Coupled Stability Model, 
Computer User's Manual" (Ref. 3). 
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SECTION II 


MODEL EVALUATION AND FORMULATION 

FEED SYSTEM HYDRODYNAMICS 
Introduction 

Liquid propellant rocket engines require a feed system to 
carry propellants from storage vessels to the combustion 
zone. Careful consideration of the feed system is required 
to produce a complete engine system capable of high perform- 
ance combined with stability. High performance systems are 
usually achieved by minimizing feed system losses, thus maxi- 
mizing the overall thrust/weight ratio. Losses, however, such 
as orifices or high pressure drop injectors provide one of 
the most direct methods of providing dynamic stability in the 
lowe r frequency range. Often then, there must be a tradeoff 
between the static and dynamic performance of the system. 
Occasionally, a feed system may be tuned to force a stable 
coupli..;) as in the use of quarterwave tubes, Helmholtz reson- 
ators or Quinke tubes. These are passive systems introduced 
to provide a resonance out-of-phase with an otherwise unstable 
system resonance. Analytical methods are helpful in (1) pre- 
dicting the dynamics of the coupled feed system, (2) providing 
a method for understanding test data, and (3) providing a 
"log^ral" test facility where, after correlation with test 
da the effect of system changes may be evaluated. 

The dynamics of a propellant feed system are concerned with 
either the determination of the feed system pressures and 
flowrate as a function of time or of characterizing the system 
frequency response. Evaluating the dynamics of a feed system 

Preceding page blank not filmed 
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requires an extension of steady-state pressure drop-flow calcula- 
tions tc include inertance and capacitive characteristics of the 
flowing system. The inertance term, as is implivd, is the ten- 
dency of the fluid to resist chanqes in its velocity due to pres- 
sure forces. Similarly, the capacitance of the fluid is the 
tendency of the fluid to resist changes in its pressure, despite 
changes in flowrate. Both the inertance and capacitance effects 
are time dependent and together describe the ability of a given 
fluid system to exhibit preferred or characteristic frequencies. 

Analytical Approach 

The feed system analysis was initially directed toward determining 
which modeling technique could best provide the data in a form 
consistent with the input requirements for the combustion dynamics 
model. As described in a later section, the appropriate output 
from the feed system model should be gain and phase as a function 
of the chamber pressure perturbation. 

Originally three approaches were considered: modal descriptions, 
lumped parameter and continuous parameter representations. It be- 
came apparent that the advantages of the modal descriptions lie in 
*heir adaptation co analog computer simulations. The lumped param- 
eter approach is limited in obtaining high frequency response for 
long line lengths (many lumps are required). However, it remains a 
plausible way for handling discrete feed system elements and can be 
tied into the continuous parameter representation which has been 
selected. 

Two methods of analyzing the continuous parameter representation 
were subsequently studied. Both yield output in the form of feed 
system gain and phase as a function of frequency. These were a 
Fourier transform method and a linear frequency response method 
(Fig. 2). By the first method, a set of nonlinear dynamic equati • 
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Flqure 2. Comparison of Feed System Modeling Techniques 
Figure c Us1 J^ a continuous Parameter Representation 








is used to generate a nonlinear transient solution. A Fourier 
transform computer program produces Fourier transforms of the 
transient data and the input signal. The ratio of these trans- 
forms, which is the system transfer function, is then evaluated 
over a frequency range and results in the system frequency 
response. Preliminary analysis demonstratpd that the Fourier 
transform program output was accurate provided that the input 
data was recorded with sufficient accuracy and that at least 
five samples were taken per cycle at the highest frequency of 
interest. That is, if 1000 Hz is the upper frequency limit, the 
time transient would have to be sampled every 0.0002 seconds. 

Employing the second method, the feed system dynamic equations are 
linearized and subsequently arranged in matrix form. The coeffi- 
cient matrix and input matrix serve as input data to a frequency 
response computer program, which then yields the required system 
frequency response. 

To determine the applicability of the two methods of dynamic charac- 
terization depicted in Fig. 2, a wave equation description of a 
four-segment baseline feed system was developed. Details of the 
analysis of the Fourier transform and linear frequency response 
methods for this baseline feed system have been given in Ref. 2. 

The results indicated that either method was capable of adequately 
describing the feed system dynamics. However, the linearized fre- 
quency response method was preferred because it provided greater 
numerical stability, was mere flexible and reauired less computer 
time. It was therefore decided that the latter approach would be 
used in the formulation of the generalized OME feed system model. 
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Development of the W at . Thanwer Equations 

Consider the differential control volume of a fluid element In a duct 
shown in Fig. 3. 



Fig. 3. Differential Pressures Developed Across the 
Incremental l^noth of a Fluid Element 


Fluid compressibility and Newton's second law leads to the following 
pair of differential equations: 


»E= _b iv = _ 2 9V 

9t P 9X 9X 


9p . 9V 6 9v 

9x p at T? at 


( 1 ) 
( 2 ) 


where 

2 2 

p = fluid pressure, N/m (lb/in. ) 

v = fluid velocity, m/sec (in. /sec) 

8 = fluid bulk modulus, N/ni^ (lb/in. 

3 3 

p = fluid density, kg/m (lb/ in. ) 

1 / 2 

c = acoustic velocity = (2/p) 

There are several ways in which to solve these equations. The solution 
method presented here follows that of Ezekiel (Ref. 4). The general 
form of the solution that satisfies either of equations ( 1) and ( 2) is 

P = Ft (t + f) + f 2 (t - |) ( 3) 

where Fj and F^ are arbitrary functions. 
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Taking the partial derivative of p with respect to x and t separately 
and substituting the results In equations ( 1) and ( 2) gives: 


-f] 

( 4) 

-e>] 

( 5) 


where 


F‘ (0 - . 


The expression for v is obtained from either equation ( .4) or ( 5): 


« ■ - F, <t + f) ♦ F ? (t - f) 


where 


z = (pB) * • 


( 6 ) 

( 7) 


Letting the subscript o denote x=0, the upstream position, and the subscript 
L denote x=L, the downstream position, and defining t = L/c as the signal 
propagation time between the two positions, equations ( 3) and ( 6) become 


P 0 * F,(t) + F 2 (t) 

P L s F 1 (t + t) + F 2 (t - t) 
zv 0 = - F 1 (t) + F 2 (t) 

ZV L = " F l^ t + T ^ + Fgft-t) 

Combining Eqs. ( 8) and ( 10), and Eqs. ( 9) and ( 11) separately, yields 
four additional relations: 


P 0 + 2v o = 2 F 2 {t) 


( 8 ) 

( 9) 

( 10 ) 
( ID 


P 0 ' 2v o * 2 F l (t) 


( 12 ) 
( 13) 
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P L + zv L = 2 ^(t-x) 


( 14) 


P L - zv L = 2 F^t+x). 


( 15) 


Eliminating the functions F 1 and F ? qives the final result as: 

L P o +ZV o] (t . T) 

[”«• - ” L ] (t . t) * >0 * «o 

Consider now Fig. 4, which depicts a generalized line segment forming a 
portion of a feed system with many such segments. 


f — 

n 

w 

n 

Flow 

Figure 4. Generalized Line Segment 
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The equations which describe the pressure and flows as functions of time 
and of each other for the generalized line segment are obtained from 
Eqs. ( 16 ) and (17 ) : 


( 18) 
( 19) 
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The expression, t - t^, indicates values at seconds before, 
and 


p n+1 ” p n+l 


+ Rjwjw 
n n‘ n 


( 20) 


w = p A v 
n M n n n 


( 21) 


Equations ( 18) and ( 19) are solutions of the wave equation, and equation 
( 20) is the flow through a nonlinear fluid resistance. Letting 




these equations can be coined to give; 

P n - a n "n ' lpn+1 * R J'V|l < 'Vl ' “nlJ (t _ T , 

p n * Vl^nKi * Vi “n * LVl * Vi VlL , 

u ’ T n-r 


(23) 

(24) 


Eliminating p R and rearranging into quadratic form results In 


Vi * ( Vl + a n> w n ' [vi * Vi VlJ 
* [V, + Vvil'Vl ' “n>] 


(t 'Vl> 


(25) 


= 0 


'“•Vi> 


which can be solved for the appropriate solution using the quadratic 
formula. The tank end parameters are obtained using a solution of 
Eq. (23 ) only. The Injector end solution is obtained using the 
quadratic formula for equation (25). 


The linear model incorporated In the Hydrodynamics subprogram utilizes 
the same basic equations, ( 23) and ( 24), but in the following linearized 
form: 
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( 26 ) 


<«P„1 - «■, <« ■ C ({ p n+l ) ' % (# Vl 0( t . ) 

<« p n+l > * “n 14 Vl> = L ({ p n> * % (s VJ (t _ T , • 


where 

(«»;*,) • (« p n +l 5 t2R n“n + l <28> 

These equations are then combined, resulting in 

«„(«"„) -< tp „> + [l ,p ntl> ,| V”n ll ‘Vl>] (t .g = 0 U9) 

<*„ + > * <• p n*l> ' [ (t p n> * «n , = ° ’ (3 °’ 


wher^ 


R = 


2 K 



(31) 


At the tank end, the term 6 p n is zero for constant tank pressure. At the 
injector end, « P n+1 is the independent variable. 


R-9807/II-9 



Injector Dynamics 

The Injector dynamics are Included by treating the 1>.„ .ct.r as a lumped 
compressible volume as shown in the figure below. 



Figure 5. Schematic of the Injector 
as a Lumped Compressible 
Volume 


The pressure in the injector manifold, pj, is related to the entering 
flow, w^, from the upstream pipe segment and the Injector flow, Wj, 
as follows: 

dp y c» 9 * 

W = VJ"g ( "n " wj) (32) 


where Vj is the injector volume and c T is the fluid sonic velocity. 


The Injector flow is controlled by the differential pressure across the Injector 
well as by the resistance and inertia of the Injector orifices. Thus, 


p i 


p c ■ R ! “I + 


1 d : 

*g 3t W I* 


(33) 


where p £ Is the thrust chamber pressure, Rj Is the injector hydr ^lic 
resistance and fc/Ag Is the equivalent inertance of all the Injector 
orifices combined, i.e.. 



n 


I 


1 

V A 1 


(34) 
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In the preceding equation, and A., are the length and area, respectively, 
of an individual Injector orifice. 


An additional factor which can have a significant effect on the response 
of the feed system to chamber pressure oscillations is injector face 
flexibility. This effect can be expressed as a change in injector voluirc 
proportional to a change in injector pressure drop: 


d V; /“ P, d P \ 
~W K [ “eft - ■ dt ) 


(35) 


Also, 



which can be rewritten as 



Combining Eqs. ( 35) and ( 37 ) gives 


d pj 

dt 


2 

* • Ct p t 

(w n ' W I> * “V7g^ 


dp I dp c 

<RT " W 


which can be rewritten as 



dp I _ C I ^ A C I P I K dp c 

ar ' V^g w n W I } + 


(39) 


This expression reduces to Eq. (32 ) if no injector flexibility exists 
(K = 0). 
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Two-Phase Flow Acoustic Velocity 

In the water-hammer equations the acoustic velocity, c, of the fluid 
appears in two places; (1) directly in the constant relating flow to 
pressure, and (2) indirectly in the time delay value, t, wh.ch equals 
fc/c seconds, where *■ is the pipe segment length. The acoustic velo- 
city of a fluid is a property of that fluid. However, its effective 
value can be reduced by the elastic walls of the flow passage or by 
the entrainment of gas and vapor in the liquid (two phase flow). 

Gas in the liquid can appear from two sources. One source is direct 
entrainment from mixing of gas and liquid in the propellant tank, 
while the other can result from the evolution of dissolved gas as 
the pressure drops along the feed system. 

Given the steady-state pressure at each point in the feed system and 
data on the solubility of the pressurant gas in the propellant as a 
function of pressure and temperature, the amount of gas in the fluid 
can be determined for each feed system segment. Then, knowing the 
amount of gas in the liquid, the effective acoustic velocity of the 
mixture may be calculated. 

Assuming isentropic compression of the gas, the change in volume of the 
gas is 

<V ' *p dp , (40) 

and for the liquid 

dV fc = \ dp (41) 

V 

Defining a constant, a = n^- 

v n 


the following relation is obtained: 




a 

a 


Kp 


( 42 ) 
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The bracketed term is the compressibility of the mixture. The density of 
the mixture can be shown to be 


aP + P, 
p m ~ (1 9 + a) 


(43) 


The acoustic velocity of a liquid in an elastic pipe is 





(44) 


Using the above expressions for density and compressibility, the acoustic 


velocity, can be written as 


c = 





p t c * 


(45) 


This expression can be used to define the acoustic velocity of a feed 
system segment with two phase flow. For an all liquid system, a = 0 
and the same equation can be used. 


In the Hydrodynamics subprogram the effect of the wall compressibility 
c f 

*erm, rr» ° n The fluid acoustic velocity is handled automatically 
te c f 

(assuming input values of ^ are provided for each feed system segment). 

However, the program does not compute the effects of two-chase flow. If 

such flow occurs in the feed system being modeled, an effective fluid 

acoustic velocity must be pre-calculated for each affected segment. 

, Cf 

Equation (45 ) above, with the ^ term set equal to zero can be used 
for this calculation. 


Simulation of Branch Lines 

In the Hydrodynamics subprogram, branched lines are handled by assuming that 
each branch has zero internal volume and that its flows are incompressib*:?. 
Thus, the pressures at the end of all segments which meet at a branch are 
set equal. The continuity of flew is then used to provide the additional 
equations in combination with the waterhanmer equations to solve for the 
overall feed system dynamic response. 
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COMBUSTION DYNAMICS 

The objective of this evaluation was to select an analytical technioue 
to describe the low and intermediate frequency dynamics of the physico- 
chemical processes leading to the combustion of storable propellants. 

The selected technique must not only be an accurate description of the 
combustion dynamics, but must be in a form that mathematically couples 
with the analytical descriptions of the feed system hydrodynamics and 
the thrust chamber acoustics. The selected technique must also include 
the effects of injection geometry, propellant properties and engine 
operating conditions. No single technique that accomplished these re- 
quirements was readily available. Instead, a combination of existing 
combustion response models, steady-state combustion models, and empirical 
correlations were integrated into a mathematical combustion dynamic des- 
cription which satisfied the technical requirements of the engireerino 
model. Before elucidating the details of the combustion dynamics formula- 
tion, a general discussion of propellant combustion is Dresented. 

Qualitative Understanding of Combustion Processes 

Propellant combustion is usually recognized as being vaporization-rate 
limited ?"d so is distributed spatially throughout the combustion chamber. 

For clarity and convenience, the combustor may be divided into a series of 
discrete zones as shown in Fig. 6 for a typical configuration. Certainly, 
transition from one zone to the next cannot be sharply defined, but is 
gradual. The positions and abruptness of these transitions are influenced 
by design variable, propel'ant combination a°d operating conditions. 

Iir.nediately adjacent to the propellant injector, the injection/atomization 
zone is least amenable to analytical description. With liquid injection 
concentrated at discrete sites, large gradients exist in all dimensions with 
respect to propellant mass fluxes and concentration, degree of atomization 
and spray dispersion, and properties of the gaseous medium. Spray droplets 
here are usually coid so the vaporization and burnina rates are low. Gases 
in this zone a^ primarily either paseous-injectants or recirculated com- 
bustion gases from the next zone downstream. 
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Figure 6. Subdivision of Combustion Chamber Into Zones for 
Analysis - Steady-State Operation 
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The primary atomization process is gradual to some extent and requires a 
finite zone length for completion, which is typically on che order of 1/2 
to 2 inches. Spray formation and its dispersion from the (approximately 
point-source) injection sites proceeds simultaneously. Frequently straight 
line ray dispersion may be a good approximation, although interactions between 
sprays from neighboring injection sites may turn the sprays. 

Completion of primary atomization and convective heating of spray droplets 
enhance vaporization rates leading to comparatively high chemical reaction 
rates in the rapid combustion zone. Upon burning, the volume of a liquid 
propellant element is increased 100-fold or more. This expansion forces 
transverse flows from high-burning-rate sites to low-burning rate sites as 
well as producing an axial acceleration. This provides some mixing but the 
sprays follow the gases only sluggishly, so spray mass flux gradients are 
primarily defined by injector- imposed geometric dispersion and interspray 
mixing. Lateral gas flows will be generated as long as appreciable spray 
flux gradients persist, but eventually they become small compared with 
axial flow velocities and the combustion field takes on a stream tube flow 
appearance. 

In the stream tube combustion zone, the flow lacks the forced transverse 
convective components that are dominant in the earlier zones. Continued 
mixing depends upon turbulent exchange between neighboring, parallel- 
flowing situations, but flow velocities here are high, residence times are 
short, and turbulent mixing is not very effective. To a good approximation, 
mixing can be entirely neglected and the two-phase flow treated formally as 
stream tubes. As sprays are accelerated and depleted, combustion rate per 
unit chamber length decays with increasing distance. Chemical reaction rates, 
on the other hand, remain high well into the exhaust nozzle. Then as the 
combustion products expand through the nozzle, diminishing pressures and 
temperatures lower the gas-phase chemical reaction rates. Two-dimensional 
flow effects are also important in the transonic and supersonic flow zones. 

For most high-combustion-effic’ency rockets, spray combustion effects are 
negligible compared with gas dynamic effects in these downstream zones. 
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Combustion Analysis 


In the past, the combustion response has been modeled with a simple time 
delay(s) (Ref. 5 through 11). This time delay represents the time required 
for the propellants to travel at their injected velocity from the point 
where they are injected to another point where they burn, and implies the 
burning is concentrated at a fixed plane some arbitrary distance from the 
injector face. The procedure outlined above is obviously an oversimplifi- 
cation of the burning process which is distributed in some fashion through- 
out the combustion chamber. 


Steady-state combustion models (Ref. 12 and 13 for example) provide insight 
to determine the droplet burning distribution as well as additional informa- 
tion required to relate the distribution to a combustion response as a func- 
tion of freauency. Combustion models are designed to march incrementally 
down the combustion chamber from a set of specified initial conditions. In 
so doing, the model calculates the rate at which the oropellants are consumed 
as a function of the axial position in the combustion chamber (burning rate 
profile). 

The analytical technique selected to describe the combustion dynamics is 
based on employing the mathematical expressions used in the steady-state 
combustion models (in particular the JANNAF DER program. Ref. 13). These 
mathematical expressions are expanded into time average and oscillatory 
components and are described in the following sections. For more detail, 
the reader is referred to Appendix B. 


Atomization Process . A very essential part of the combustion fi 'Id initial- 
ization is the assignment of propellant spray droplet sizes and flowrates. 
Analytical descriptions of the atomization process are not available but 
empirical correlations that relate droplet diameter to injector geometry 
and flow conditions are available (Refs. 14, 15, and 16). For like-doublets, 
one empirical relationship is (Ref. 14). 


D d - 4.85 x 10 4 Vj- 0 - 75 (P c /Pj)-° W dj 


0.57 


(46) 
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where Vj Is the liquid jet velocity and dj is the liquid jet diameter 
at the atomization plane. (For steady-state analysis, the velocity is 
the injection velocity and the diameter is the orifice diameter.) 


For purposes of the current analysis, the atomization process is des 
cribed by 


K(dJ a 
J x=x. 


<*/ 

imp J x **1mp 


( 47 ) 


where x^ is the location of the atomization plane or the impingement 
point. Expanding Eq. 47 into tira-averaged and oscillatory parts, yields 
the oscillatory droplet die.iieter 



In order to evaluate the oscillatory droplet diameter, che oscillatory 
liquid jet diameter and velocity (and therefore the jet flowrate) are 
required at the atomization plane. Therefore, the dynamics of the fluid 
from the injector to the atomization plane, the Fenwick and Bugler Klystron 
effect (Ref. 17) is required. The dynamics of the propellant transport 
process to any location in the chamber are described by the continuity 
and momentum equations (for constant density): 


_ 3 _ 

3t 


<v 


& (A J v ■ 0 


(49) 


3 

at 


(A j V * 


_ 9 _ 

ax 


(a j V> 


0 


(50) 


Expanding the preceding equations into time average and oscillatory parts 
and integrating between the injector face and any location in the chamber 
yields 
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(51) 


v. 


J. = J“*/v, 




(52) 


(53) 




where u> is the complex frequency and the oscillatory injection rate, (m,) 

J i n j 

is determined by the feed system dynamics. Equation 53 Is the oscillatory 
flowrate at x and is usually referred to as the Klystron effect (Ref. 17). 
The Klystron time delay, t^, is therefore given by 


c = 


(54) 


Considerable amplification of the injector face flow oscillations are 
possible when the Klystron effect is present and could explain the periodic 
burst of acoustic resonances called resurqing and the steep-fronted waves 
seen in low and intermediate frequency instabilities. 

Droplet Vaporization . Theories of droplet combustion (Refs. 12, 18, 19) 

are available which may be used to evaluate the extent of coupling between 
droplet burning rate and local pressure and velocity luctuations. In 
general, droplet burning is enhanced by increased turbulence levels or by 
periodic directional variations in velocity, because droplets are relatively 
heavy and resist following gas streamlines. 

Calculation of the spray heating and vaporization is usually accomplished 
through specification of the corresponding individual droplet processes 
and sunmation over all the droplets that constitute the spray(s) being 
analyzed. The calculation of single droplet evaporation is usually based 
on a spherically symmetric model of simultaneous heat transfer and mass 
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transfer across the gas side boundary, o** film, separating the liquid 
droplet from the surrounding hot contoustion gas. Forced convection and 
resultant nonspherical transfer processes are accounted for through 
empirical Nusselt number correlations for both heat and mass transfer. 


Based on the derivations presented in Ref. 12, the droplet vaporization 
rate for the k— droplet group size may be expressed as 


m. 


( 6 ) Z k k f. Nu H. m k 
k k 

VaPk " ^ C d 2 'n 
l k d k C v. 


(55) 


where c Nil P MW D 
p v k m k v k v k 

Z k 5 H ' Nu u ft T. 
f k H k u f k 


in 


VP - Pk 


(56) 


The fuel or oxidizer droplet spray continuity equation is 

m. 


I 37 (A c k V = - A I p 


van, 


k m, 


(57) 


where is the spray density, v k is the droplet velocity and m^ is the 
mass of a single droplet. 


For steady-state combustion models, the preceding equations are numerically 
integrated allowing the droplet diameter, , to vary along the length of 
the combustor. In order to simplify the integration for stability analysis, 
the droplet diameter was held constant and the droplet number flowrate was 
assumed to vary (for mere detail see Appendix B ). Combs (Ref. 20) has shown 
that changing from a variable droplet diameter to a variable droplet number 
flowrate yields approximately the same results for steady-state vaporization. 
Therefore, substituting Eq.C5 into Eq. 57 and lettinq 
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(58) 



Integrating Eq. 62 between the start of vaporization (x^) and any general 
location (x) and substituting the resulting expression into Eq. 62 yields 
the spray vaporization rate 
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(64) 


(65) 



where 


(v s ) 


x=x,, 


s ~ 




x-x. 


T 


s 


T 

s 



( 66 ) 


(67) 


and the oscillatory spray droplet diameter (0 $ ) is qiven by Ea. 48 and 
the oscillatory flowrate is given by Eq. 53. The above formulation 
results in a linear oscillatory vaporization model similar to Crocco's 
n-T model (Ref. 6). The formulation includes the effects of: (1) dis- 

tributed energy release, (2) oscillations in the injection rate, (3) oscil 
lations in droplet diameter, (4) oscillations in droplet temperature, 

(5) gas pressure and velocity oscillations, and (6) oscillations in the 
local mixture ratio. 


Nuss el t Number . It may be observed that one of the dominant terms in both 
the expressions for the average and oscillatory time delay is the Nusselt 
number. The Nusselt number, tor longitudinal modes, is (Ref. 21). 


Nu 


H 


2.0 


+ 0.6 Pr 


1/3 




(68) 


In order to evaluate the oscillatory Nusselt number, the oscillatory droplet 
spray velocity is reouired. The droplet spray velocity can be obtained from 
the drag equation. 


m 


dv s 

ar 



Iv-v s l(v-v s ) c D 


(69) 
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Let 



where 


a 

S 



( 70 ) 


(71) 



and the oscillatory Nusselt numbor can be written as (for more detail, 
see Appendix B) 



Calculations have been made which indicate that, for large droplet diameters, 
the average and oscillatory Nusselt numbers are quite sensitive to pressure 
and velocity oscillations. Therefore, the Nusselt number can have a signi- 
ficant effect on engine stability. 
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Droplet Heat Transfer Bloc x age Term. The oscillatory combustion t^e 
delay given by Eq. 67 requires the evaluation of the he*t transfer 
blockage term (Z $ ) which is related to the ccmbustion gas and liquid 
vapor properties by Eq. 56 . Because the vapor pressure (P v ) at the 
droplet surface is related to the droplet temperatur tne blockage term 
also depends on t » oscillatory droplet surface temperature inside the 
droplet which is given by: 


J_ (p c t ) * -l — 
9t ' f v^ V 9r 


r 2 k 


3 A 

eff » 8r 


(76) 


Therefore, the oscillatory heat transfer rate to the droplet can De 
related to the oscillatory droplet surface temperature by 


Q 


s 



(77) 


where Ry is given in Appendix B. The droplet heating rate can also be 
written is (Ref . 12) 


Q s ■ Z $ kf % 


A H, 


< T - V . _ 

(e Zs - 1) X 


vap $ 


(^ s ) 


(78) 


Assuming that 



= 0 


(droplet at "wet bulb" temperature) 


( 79) 


and 


= = R p 

n 


! P 

it/ 


(80) 
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etc. for the other variables, the response factor for the heat transfer 
blockage term can be related to droplet and gas properties and flow 
conditions (see Appendix B ). 

Examination of the response factor for the heat transfer blockage term 
indicates that this term can have a dominant effect on stability if the 
frequency is high and the effective thermal conductivity is low. 

Generalized Vaporization Rate Expression 

In order to maintain generality in representing the combustor dynamics, 
the spray vaporization rates (fuel and oxidizer) can be written as: 



From the preceding sections, and Appendix B , the combustion coefficients 
(c.| through c^g) can be defined in terms of injector geometry, flow condi- 
tions, etc. 
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CHAMBER DYNAMICS 

During certain periods of a rocket engine's operation, conditions within 
the combustion chamber are time variant, i.e., the operation is not steady 
with respect to time. In the following sections, prime interest is focused 
on abnormal transient operation during unstable combustion, i.e., pressure 
oscillations in a combustion device which are sustained by the combustion 
process. Start and stop transients are not considered. 

Classifications of Instability 

The deviations from steady-state combustion which occur during unstable 
burning depend upon the kind of instability experienced. Liquid rocket 
instabilities are classified according to their dominant time-varying 
processes. They may be divided initially into two categories, depending 
upon whether the instability oscillation wave length is long or short 
compared with the chamber dimensions. 

If the instability wave length is considerably longer than the chamber 
length and diameter, pressure disturbances propagate rapidly through the 
combustion space compared with rates of change due to the instability. 

As a result, wave motion in the chamber may be neqlected and chamber pres- 
sure can be considered to vary only with time but not to vary spatially 
(i.e., P c is a lumped parameter). These instabilities depend upon a fluid 
mechanical coupling between the propellant feed system(s) dynamics (fluctuat- 
ing injection rates), the propellant combustion rates (delay times), and the 
combustion gas exhaust rates (pressure relaxation). Such instabilities 
can be further subdivided into various categories depending on the extent 
of wave motion in the feed system. 

The breakpoint at which chamber wave motion becomes important is not abrupt. 
In reality, chamber wave motion is always present and, in effect, lumped 
chamber instabilities a^e really "zero order mode" limits of more general 
wave motion instabilities. In practice, it is found that the chamber gases 
can be considered to act as a lump until the frequency of osci 1'ation exceeds 
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roughly one-fourth of the frequency of the lowest chamber acoustic 
resonance mode. At and above such frequencies wave motion becomes 
important and cannot be neglected in analysis. Chamber wave motion 
instabilities are characterized by wave-length of the oscillatory 
motion being comparable to the chamber dimensions. As with lumped 
chamber instabilities, the driving energy comes from oscillatory spray 
combustion. With wave motion instabilities, however, in addition to 
the effects of injection rate fluctuations, there is the combustion 
response of burning propellant sprays as they are disturbed by passage 
of a pressure wave through them. Wave motion may increase local burning 
rates by any of several mechanisms: (1) a pressure effect on the drop 
vapor gas phase burning rates; (2) enhanced mixing between gases and 
between sprays and gases; and (3) increased spray gasification rates. 
Increased spray gasification may be due to transient increases in con- 
vective flow velocities, to increased temperature or concentration 
gradients, and/or spray droplet shattering. The instability amplitude 
depends upon the magnitude of the resDonse, and vice versa; typically, the 
interacting processes are driven to a limit represented by abrupt, es- 
sentially complete consumption of the propellant sprays. This direct 
response can be so great that injection rate fluctuations may be of 
secondary importance. As a result this class of instability can also 
be further subaivided as to the importance of feed system coupling. In 
the absence of feed syster., coupling, the instability is referred to as 
"classical acoustic instability." Only feed system coupled instabilities 
are considered in this program. 

Analytical Approach 

Two methods of approach were considered for solving the chamber dynamics. 

"he f : rst method used a linear lump chamber coefficient. This method is 
valid only at low frequencies (less than 500 Hz) and results in a set of 
nonlinear algebraic equations to be solved. 

The second method employed a first-order perturbation model to define the 
chamber frequency and growth coefficient along with the oscillatory pressure 
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distribution in the chamber. This method is valid for all frequencies 
of interest in the present program (10 to 1000 Hz). For the oscillatory 
variables, solutions of the form $ = <j>' e” 1ljt , where u> is the complex 
frequency, were assumed. These forms yielded a set of nonlinear dif- 
ferential equations which were numerically integrated between the injector 
face and the nozzle inlet plane. Using iteration techniques and the 
requisite boundary conditions at the injector and nozzle inlet plane, 
the chamber frequency and growth coefficient were obtained. 

Consideration of the degree of complexity in solvinq the governing eoua- 
tions by each of the above methods, as well as the range of validity of 
each approach* resulted in choosing the first-order perturbation models as 
the best method for describing the chamber dynamics. In the following 
paragraphs, the solutions to the first-order perturbation model stability 
equations are p* ~d without showing their derivations. The reader 

is referred to x C for a complete derivation of the chamber model 

equations. 

First-Order Perturbation Model 

The perturbation technique is a useful mathematical tool to simplify non- 
linear partial differential equations. Assuming a variable <*> can be ex- 
pressed as: 

$ = <j> + $ (82) 

_ ** 

where 4> and $ represent the time-averaged and oscillatory components, 
respectively, and letting $n << $n, results in a set of equations of 
reduced complexity. 

The assumptions used in the derivations of the chamber model equations are: 
(1) ideal gas flow is a valid state equation, (?) dilute sprays occupy a negli- 
gible fraction of the chamber volume; (3) the spray can be represented 
by a finite number of dropsize groups; each dropsize group contains a 
large number of locally identical drops; and, each size group constitutes 
a separate liquid phase and exchange terms between liquid phases are not 
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Included; (4) drag contributes only kinetic energy to the spray energy 
equation; (5) secondary "shear" breakup of drops is not included; 

(6) no body forces; (7) onc-dimer c ional axial flew; (8) diffusion, 
thermal and viscous gradients are negligible; and (9) droplet drag 
forces and heat transfer to the droplets are negligible. 

Based upon these assumptions and following the perturbation technique 
described above, the equations for longitudinal instabilities can be 
written as (see Appendix C). 


(a) Gas Continuity 
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(b) Gas Momentum 


v .(dii) +vl dv + v__ dvT_ p' dp 4 _E ' ijE + — E gEl 

' ' c* dx Cj? dx — „ 2 3x — \ 2 Bx - ^ 2 dx 


p c 


0 


p c 


0 


P c. 


(c) Equation of State 
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(e) Gas Mixture Ratio 
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The preceding set of ordinary differential eauations are numerically 
integrated once the complex frequency is specified, between the injec- 
tor face and the nozzle inlet plane. The method of calculating the 
complex frequency for the perturbation model, based on nozzle admit- 
tances calculated from upstream and downstream variables, is discussed 
in the Computer Model Development Section under the Engineering Model. 
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SECTION III 

MODEL DEVELOPMENT AND SOLUTION 


FEED SYSTEM HYDRODYNAMICS 
Introduction 

Prior to the development of the generalized hydrodynamics model , a trip 
toMcDonnell Douglas/St. Louis was made to obtain design and operating 
mode data for the OMS, PBK and RCS feed systems. The data included 

layouts and dimensions for all the lines in the propellant flow paths, 

as well as information related to propellant properties and system operat- 
ing pressures, for the various system operating modes. 

The data were found to be satisfactory for definition of the generalized 
feed system model such that suitable representation of the OMS feed system 
could be employed. 

Generalized Feed System Model 

Based upon the design and operating mode data for the OMS, PBK and RCS 

systems, a generalized feed system model was developed f*r use in the 

Hydrodynamics subprogram and is shown in Fig. 7. The system is comprised 
of 30 individual line segments, each denoted in Fig. 7, as the lines 
between the black dots. As described in Section II, a continuous parameter 
representation of each line segment is obtained through the use of separate 
sets of waterharrmer equations. Each line segment can have a different line 
length, area, wall compliance, fluid acoustic velocity and resistance, and 
hence can model a wide variety of feed system components by merely choosing 
the appropriate values from these parameters. Also included in the generalized 
model are lumped parameter descriptions of two injectors (designated 0 and F 
on Fig. 7). Parameters for the injectors are volume, resistance, inertance, 
fluid acoustic velocity and face flexibility. A list of feed system compo- 
nents which may be included in the generalized OME feed system analysis are 
given in Table 1 . 
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Figure 7. Generalized OME Feed System Schematic 
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TABLE 1 . FEED SYSTEM FEATURES MODELED BY THE 
GENERALIZED OME HYDRODYNAMICS MODEL 


• Propellant Tanks 

• Propellant Feedlines 

• Gas/Propellant Heat Exchanqer 

• Orifices 

• Screen Filters 

• Flow Control Valve (variable area cavitating venturi) 

• Propellant Shutoff Valve 

• Flexible Bellows 

• Propellant Injector 

• Rocket Engine Thrust Chamber 

• Ablative Cooled 

• Regeneratively/Film Cooled 

• Rocket Engine Nozzle 

• Propellant Accumulators 

• Propellant Acquisition System 
o Propellant Properties 

• Variable Density and Temperature 

• Pressure 

• Quality 

• Pressurizing Gas Saturation 

• Structural Incuts 
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The system of 57 equations describing the generalized feed system de- 
picted in Fig. 7 is provided in Appendix A. It should be noted that 
the Hydrodynamics subprogram solves the complete system of equations 
each time it is called. Thus the frequency response of the entire 
system is calculated each time. It has been shown, however, that simpler 
feed systems, representing only a portion of the Fie. 7 schematic, can 
be modeled by merely assigning values to the parameters of the unneeded 
line segments which will exclude them from having any effect on the system 
frequency response (Ref. 2). This is accomplished automatically by the 
Hydrodynamics subprogram via the assignment of very large resistances and 
Vfe*'y short lengths to all line segments for which no data is entered. 

Solution Method 

The system of equations describing the generalized feed system is solved 
using standard matrix techniques. A subroutine (FRESP) in the Hydrodynamics 
subprogram (HYDRDY) solves for the variables X. in the following relationship 

[C] • (X) = a . LY J 

where {X} represents the column matrix of independent variables (pressures 
and flowrates), and Y is a single input variable that represents a unit 
value of the injector end combustion chamber pressures. The matrix [a] 
then relates the specific pressure input to each applicable equation that 
contains combustion chamber pressure ( [a] may contain both static and 
dynamic terms.) The matrix [c] is simply the coefficients of the linear 
differential equations that represent the physical system. The values of 
the coefficients for the [a] and [CJ matrices are computed by the sub- 
routine HYDRDY. 


The FRESP matrices can be expressed as: 


C ijk S 


k-1 


{ y - 


a 1k S 


k 1 


( 88 ) 
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with the differential operator defined as S = Jw, where J = and 
u is the frequency. The matrices may be broken down to provide real 
matrices and imaginary matrices. 
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The matrices are multiplied and then solved for [xj using a standard 
Gaussian elimination procedure for solving linear equations. The [x^] 
solution is still separated into real and imaginary components, and are 
simply combined to form a vector for each variable. The procedure is 
repeated for each frequency being considered. 

Structure of the Hydrodynamics Subprogram 

A functional block diagram of the Hydrodynamics subprogram (HYDRDY) is 
shown in Fig. 8. HYDRI"' is called by the main program to calculate 
the frequency response characteristics of the feed system. Its functions 
include (1) reading of input data describing the physical attributes of 
the feed system components, (2) generation o. a matrix of linear dif- 
ferential equations representing the complete feed system, (3) solution 
of the feed system equations to yield the amplitude and phase response 
of all feed system pressures and flowrates as a function of chamber pres- 
sure oscillation amplitude and frequency, and (4) generation of tabulated 
output of injector flowrate frequency response for use by the main program. 

A basic assumption of subroutine HYDRDY is that the feed system being 
modeled can be represented by the generalized schematic of Fig. 7 (or 
by some portion of this hematic). This assumption is necessary be- 
cause HYDRDY sets up and solves the complete set of simultaneous eouations 
representing the Fig. 7 schematic. By assigning very hi qh resistance 
and very short length values to any of the 30 numbered line segments of 
the generalized schematic, those segments can effectively be excluded from 
having any effect on the frequency response characteristics of the rest 
of the system. With this approach a wide variety of feed systems can be 
modeled with no changes to the program other than the input data. 

Input control variable IR directs the reading cf subroutine HYDRDY input 
data. If IR is zero or less, the program assumes that all required data 
has previously been read and the data read function is bypassed. If IR 
equals one or greater, provision is made to read in the appropriate data 
for the line segments and injector(s) to be included in the analysis. 
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Figure 8* Block Diagram of Hydrodynamics Subprogram 





Control variable INPHYD determines the procedure by which the oscillatory 
injection flowrates arc calculated. If INPHYD is less than or equal to 
one, the hydrodynamic calculations are performed for one iteration step or 
frequency. If INPHYD is qreater than one, the hydrodynamic calculations 
are performed for a range of frequencies and the results stored on taoe. 

Simultaneous solution of the 57 linear waterhammer and injector equations 
describing the complete Fig. 7 generalized feed system, at each specified 
input frequen.y, yields the oscillatory amplitude and phase response of 
all pressures and flowrates in the feed system to inputs via chamber pres- 
sure oscillations at that frequency. It should be noted that although out- 
put from a single call to HYDRDY contains values for both "oxidizer" and 
"fuel" oscillatory injection flowrates (at one or more frequencies), the 
output values actually refer to the "0" and "F" injectors of the Fig. 7 
schematic. Thus, unless both oxidizer and fuel feed systems can simul - 
taneously be modeled With the Fig. 7 layout, it is necessary call 
HYDRDY twice - once for the oxidizer feed system and once for the fuel 
feed system. 

A complete description of the Hydrodynamics subprogram, as ..'ell as 
detailed information related to all input variables, is qiven in the 
CME/Feed System Coupled Stability Model, Computer User's Manual (Ret. 3). 
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COMBUSTION DYNAMICS 


Str ~ture of the Combustion Dynamics Subprogra m 

A block diagram of auxiliary programs required to evaluate particular 
combustion parameters is presented in Fig. 9. Equilibrium properties 
of the combustion gases are evaluated using the NASA ODE computer 
prog- am, or an equivalent program. An existing Distributed Energy 
Release comoustion model (JANNAF computer program DER) which uses 
the combustion gas properties, chamber and nozzle geometry, injector 
geometry, and engine operating conditions is used to calculate 
distributed combustion parameters. The distributed vaporization 
rates calculated by the DER computer program are used to evaluate 
the average droplet vaporization "time delays", Nusselts numbers, 
and other important steady-state parameters required by the combustion 
dynamics subprogram. Output from the DER program, oscillatory 
injection flowrates, engine operating conditions, propellant properties, 
chamber geometry, and the Klystron constants are input data for the 
combustion dynamics subprogram. 

A functional block diagram of the Combustion Dynamics subprogram 
(C0MBDY) is shown in Fig. 10. C0MBDY is called by the main program 
to calculate the combustion coefficients required by the chamber 
dynamics routine. Its functions include (1) reading of input data 
describing the physical attributes of the combustion process, (2) 
solution of the combustion equation to yield combusti;n coefficients 
and (3) generation of tabulated output of combustion coefficients as 
a function of frequency for use by the chamber dynamics routine. 
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Figure 9. General Structure of Combustion Dynamics Model 
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Input control variable IR directs the reading of subroutine C0MBDY input 
data. If IR is less than or equal to zero, the program assumes that 
all required data has previously been read and the data read 
functions are bypassed. If IR is greater than or equal to one, 
provision is made to read in the appropriate data for evaluating 
the combustion effects to be included in the analysis. 

Control variable INPC0M determines the procedure by which the combustion 
coefficients are calculated. If INPC0M is less than or equal to 
one, the combustion coefficient calculations are performed for 
only one frequency. If INPC0M is greater than one, the combustion 
coefficients calculations are performed for a range of frequencies 
and the results stored on tape. 

A complete description of the combustion dynamics subprogram, as 
well as detailed information related to all input variables, is 
given in the OME/Feed System Coupled Stability Model, Computer Manual 
(Ref. 3). 

Combustion Coefficient Evaluation 


The main purpose of the combustion dynamics subprogram is the 
evaluation of the combustion coefficients which appear in the generalized 
vaporization rate expression(p. 11-28). Basically, the coefficients are 
evaluated by substituting the oscillatory equation for the combustion 
processes into the oscillatory vaporaization rate expression and then 
comparing the resulting equations with the generalized vaporization 
rate expression. 

For example, the oscillatory spray velocity is 
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Si nee 


(f) ' H £) 

' 'inj ' p 'inj 


(93) 


*he oscillatory spray velocity can be written as 


e 1 “ x k/' v s 6. , 
s inj. 


m . . 

\ P / mj 


(94) 


To eliminate the time dependency from the preceeding expression, recall 

0 = 0 ' e ~ i<,,T 


(95) 


therefore 


i “ x k/ ; s r 

e s * G fnj s e s 


(3 


1U)T 


inj 


G . 

mj. 


-ito (t- x k /v s ) 


(96) 


inj 


The bracketed term in the preceeding expression is just the injector end 
pressure evaluated at earlier tin 
spray velocity can be written as 


pressure evaluated at earlier time ( x k $ /v s ); therefore, the oscillatory 


Gl " j (f) 

x=o 


(97) 


where / P \ includes all time phases which is consistent with the method 
Wx-o 

which is used in the chamber dynamics equations. 
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Using similar methods to evaluate the other oscillatory combustion processes, 
the generalized combustion coefficients are 
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(107) 

('.08) 
(109) 

It should be noted that the oscillatory heat bloc’ age term has been 
neglected based on work presented in Refs. 22 and 23 which indicates 
that this term is not important at low frequencies. 

Combustion Characterization Test Plan 


where the subscript s denotes the fuel or oxidizer and 



One of the inherent analytical limitations of the Feed System-Coupled 
Stability Model is to be found in the determination of the complete set 
of combustion parameters required as input data to the combustion dynamics 
subprogram. The combustion dynamics model has been formulated to permit 
variables obtained from real combustion systems to be related to the mathe- 
matical model. The Klystron constant, wnich results from the analysis of 
the oscillatory fluid dynamic properties from the injector to the atomiza- 
tion plane (Ref. 17), is one example of a parameter that must be experimen- 
tally determined for the particular combustion system being investigated. 

Originally, the experimentation needed to obtain the combustion parameters 
was to be conducted by NASA concurrently with the present program's effort. 
To aid in these determinations, and as partial fulfillment of the require- 
ments of Contract NAS9- 14315, a "Combustion Characterization Test Plan" 
(Ref. 24) was written which defined in detail the test objectives, test 
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hardware, test procedures, and data requirements necessary to obtain em- 
pirical characterization of stability-related propellant injection param- 
eters and sensitive operating conditions for OME-type combustor hardware. 

The propellant injection parameters to be included were nozzle admittance, 
injection admittance. Klystron response, and total combustion time lag. A 
test series was proposed that involved four subscale injector configurations 
(three like-doublet units and a triplet injector) in conjunction with a heat 
sink combustor using the ^C^/MMH propellant combination. Details of the 
series of 36 tests were outlined for the determination of stability limits 
for the four injector configurations. A series of 48 tests was also recom- 
mended to determine the combustion response of two of the injectors over a 
range of chamber pressure, mixture ratio, and frequency. 

The above plan was subsequently changed due to delays in conducting the ex- 
perimental program. As a result, the combustion parameters required by the 
present computer model for use in model verification were estimated using 
the best available experimental data. 
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CHAMBER DYNAMICS 


Structure of the Chamber Dynamics Subprogram 

A block diagram of auxiliary programs required to evaluate the gas flow 
variables and the chamber dynamics is presented in Fig. 11. An existing 
Distributed Energy Release combustion model is used to calculate 
distributed combustion parameters which are employed in defining the 
steady-state flow variables. The combustion dynamics subprogram supplies 
combustion coefficients for evaluating the oscillatory distributed energy 
release along with other important combustion parameters, such as the 
vaporization time delays. Klystron constants, etc. Additional input 
required by the chamber dynamics subprogram are the chamber geometry and 
the nozzle admittance based on downstream conditions. 

Output from the chamber dynamics '‘.j_,.**ogram are oscillatory chamber gas 
variables as a function of dist_ :c from the injector face and a comparison 
between tne nozzle admittance calculated based upon both upstream and 
downstream conditions. A complete description of the combustion dynamics 
subprogram, and other supporting routines, is given in the OME/Feed System 
Coupled Stability Model, Computer Manual (Ref. 3). 


Solut ion ^■ ‘thod 


Th.* ordinary differential equations describing the oscillatory solution 
a so;v;**i using a second order implicit finite difference method. This 
method has the advantage of being simple to implement and modify, as 
well as being unconditionally stable for systems of equations which do not 
have exponentially growing solutions. The method as applied to the first 
order system 

y* = Ay + g (no) 

where Y and g are nxl vectors and A is an nxn matrix is as follows: 

y i+l = y i * T A i+l/2 (y i 4 y i+l ^ + 9 i+l/2 (111) 
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Here, the subscript 1 refers to the I'th mesh point in the finite difference 
scheme, e.g., Xj s x Q + iAx. The approximate the Y vector at x^. That 
is y 1 “ Y i = Y(x^). The subscript 1+1/2 refers to evaluation at x^ + ax/2; 

e *9* * A f+i/2 = A ( x i + ^ 2 ’)* 

That the above method leads to a second order approximation (error is 

O 

proportional to ax ) can be shown as follows: 

Solving for y^ +1 yields 

y *+l = (I - Y A i+i/2^ 1 (I + Y A i+1/2 ) y i + AX (I - Y A i+l/2^ 1 g i+l/2 

( 112 ) 

Without loss of generality, assume i = 0. 


From the two expansions 


Y s Y + — Y* + 

Y 1 ' 1/2 2 Y l /2 


y « y -My 1 + 
Y 0 T l/2 2 T l/2 


4 YV/2 + o(ax 3 ) 
IT Y l/2 + °^ Ax3 ^ 


(113) 

(114) 


the following are obtained 

Y^ = Yq + Ax Yy 2 + ) (115) 

and 

Y 1/2 * ( Y 0 + V' 2 + o(ax2) 

Let y Q - Y q ; it is necessary to show y-j 
second-order accuracy. 

From the differential equation (110), one obtains 

g l/2 = Y l/2 " A l/2 Y l/2 (117) 


(116) 

3 

= Y^ + o(ax ) in order to demonstrate 
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Substituting (117) into (112) and noting Yq = Yq yields 


y l = (I - f A 1/2 ) _1 (I + f A 1/2 ) Y q + Ax (I - f A 1/2 ) _1 (118) 

^ Y l/2 ' A l/2 Y l/2^ = ^ ‘ T A l/2^ 1 | Y 0 + Ax Y 1/2 + A l/2 Ax ^ 1/2 Y 0 * Y l/2^| 

Using (115) and (116) gives the result (119) 

y l = Y 1 + (I * T a i/ 2 ) ~ 1 o(ax3) {120) 

= Y 1 + o(Ax 3 ) 

Consider now the stability of the finite difference formula ( ill) for systems 
which do not have exponentially increasing solutions; that is, the real 
part of each of the eigenvalues of A is negative. To prove that they are 
stable for this situation, define the error * Y,. - y^ and consider the 
two equations given by (111) and 

l+l u 2 ft i+l/2' u 2 A i+l/2 ; T i 

ax (I - f A i+1/2 ) -1 g. +1/2 + o(ax 3 ) , (121) 

the latter resulting from (120). Subtracting (HI ) from(121) yields 

e i+l ^ " T A i+l/2^ ^ + ~Z A i+l/2^ e i + °( Ax ^ 

Let B = -y- A^ + y 2 - The method is stable if and only if the matrix (I-B)"^ 

(I+B) has a spectral radius less than one, for this would produce (Ref. 25) 


lim 

n-*o 




(123) 
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Since the eigenvalues of (I-B) -1 (I+B) are just equal to 0+B)/(1-8), 
where 6 Is an eigenvalue of B, the spectral radius of (I-B)"^ (I+B) 
is just 

max |(l+8)/(l-e)| (124) 

8 

For this to be less than one, 

1 1 +6 1 < 1 1 - 8 1 (125) 

for all 8. This implies 

1+8 + 8 + 8 ? < 1-6-? + 8? (126) 


or 


6 + 8 < - (8+ e) 

(127) 

Real (e) < 0 

(128) 


Since 6 = ^ » where a is an eigenvalue of A, the method will be stable if 
all the eigenvalues of A have real parts less than zero, that is, the 
solutions to (110) are not exponentially increasing. 
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ENGINEERING MODEL 
Structure of the Program 

The general structure of the Engineering Model is depicted in Fig. 1? . The 
decision to adopt this structure was based on a trade-off o* setup time, 
storage capabilities, and solution time. An equilibrium gas property 
program similar to NASA ODE computer program (Ref. 26), and the DER combustion 
model (Ref. 13) must be executed external to the stability program in order 
to obtain important input information. General input data to the program 
including geometric factors, engine operating conditions, propellant 
properties, starting location for vaporization, estimated complex frequency, 
and control variables are r st input to the program. The control program 
will then execute the nozzle admittance and hydrodynamics subprograms each 
time through the program or calculate the admittance and/or oscillatory 
injector flowrates as a function of f quency and store the results on tape. 
Input data required by the nozzle admittance and hydrodynamic subprograms 
are read in the first time these programs are executed. Next, the combustion 
dynamics and steady-state subprograms are called. The combustion dynamics 
subprogram inputs combustion data the first time it is executed end the 
steady-state subprogram is only executed the first time through the program. 

The chamber dynamics subprogram then numerically integrates the chamber 
dynamic equations from the injector plane to the nozzle inlet plane and 
calculates the nozzle admittance based on upstream conditions. This 
admittance is then compared with the nozzle admittance based on downstream 
conditions. Numerical methods are then used to calculate a new estimate 
for the complex frequency and the above procedure repeated until the two 
admittances agree within a prescribed limit. 

Output of the model includes a prediction of the oscillatory frequency and 
information related to the stability or damping coefficient (decrement) of 
that frequency for the particular propulsion-fed system. In addition, at 
each solution point, the amplitude and phase angle for the oscillatory pressure 
ratio, velocity ratio, temperature ratio, and mixture ratio are given as a 
function of distance from the Injector face. 
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Determination of Complex Frequency: Numerical Solution He thod 

The complex frequency, u> , is determined such that the boundary condition in 

the nozzle is satisfied. Specifically, the admittance is reauired to be 

continuous across the interface between the combustion zone and the zone 

immediately downstream of the combustion zone. In the downstream zone, 

the nozzle admittance, » is computed from a nozzle admittance program. In 

the upstream combustion zftne, the nozzle admittance, is computed from the 
oscillatory flew parameters determined by the chamber dynamics. The com- 
plex frequency must be such that 

A n , A n (129) 

N D U 


Let ai = x + iy and F 
x and y such that 



u + iv. The numerical problem is to find 


u(x,y) = 0 (130) 

v(x,y) = 0 (131) 


Several methods were considered for solving this system of equations. 
Because F is not an analytic function of u>, the complex form of the Newton- 
Raphson method may not always work. On the other hand, one could use the 
two-dimensional form of Newton-Raohson (Ref. 27)> but since the derivatives 
of u and v with resoect to x and y must ta computed numerically, the two- 
dimensional Newton-Raphson method will require three functional evaluations 
of F at each oi, i.e., (x,y), (x + Ax,y), and (x,y + Ay). Alternatively, a 
far more efficient method is to use the two-dimensional form of the secant 
method (Ref. 27) since this does not require the evaluation of any deriva- 
tives. Specifically, this method approximates the u and v surfaces with 
linear functions u L and v^ (planes) based on three previous guesses for us 
[(x-i.y-i). (x 2 ,y 2 ), (x 3 ,y 3 )]. The next guess for w, (x 4 ,y^), is determined 
from the equatio s u L (x^,y^) = v L (x^,y^) = 0. Tne new value of oj then re- 
places one of the previous three values, normally the one with the largest 
error as measured by the absolute value of F(x.,y.), and the iteration 

sJ J 
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is continued until convergence is reached. The actual equations for the 
above process take the following form. Let w. = x. + iy., u. = u(x.,y.), 
and v. = v(x.,y.) for j = 1,2,3. J J 3 J 


1. Determine it., j = 1,2,3, such that 

J 


"l + *2 + ”3 = 1 

(132) 

*1 °3 + V 2 + *3 U 1 = 0 

(133) 

"1 v 3 + ff 2 v 2 + Vl = 0 

(134) 

2. Compute uj^ = + ^2^2 + 7r 3 w ’| • 

(135) 


3. Compute u^ and v^, based on and y ^ 


4. Test for convergence, i.e., require that min |uk - u^|/|w 4 l < 

and |F 4 | < e 2 . 3 

If the process has not converged, continue with steps 5, 6, and 7. 

2 2 

5. Determine the j between 1 and 3 such that u. + v. is maximum. 

J 3 

6. Replace u>., u., and v. by to., u., and v., respectively. 

J J J “ *T • 

7. Go to 1. 


Operationally, steps 4 and 5 may be altered to replace the to's cyclically, 
i.e., uk -*• -*■ v^ ->■ In fact, the computer program as 

written alternates between these two procedures every three iterations in 
order to avoid any possible cycling that may occur. 

The above algorithm has been found to be very efficient when the first three 
guesses are relatively near an actual solution. The difficult problem was to 
develop a searching algorithm which determines the regions in the uj plane 
where lutions exist. 

2 2 

One possible procedure would be to utilize the fact that the surface u + v 
has an absolute minimum at each solution. Using any reasonable value of to 
as a first guess, one might be tempted to employ a gradient, or modified 


R-9807/ III -25 



gradient, method to march along the surface until one came near a relative 

2 2 

minimum. Unfortunately, this procedure fails because the surface u + v 
has many relative minima which are not actual solutions. The reason for 
the large number of relative minima (and maxima) for this surface is undoubt- 
edly due to the coupling between the combustion processes and the feed system 
oscillation in conjunction with the very rapid change of the feed system 
response as a function of frequency. The searching algorithm must be able to 
discriminate between those relative minima that are not solutions and those 
that are. Such a procedure was developed for this program. It takes 
advantage of the fact that a large portion of the computations required are 
only a function of the real part of w, i.e., they use x as an independent 
variable and do not depend upon y. Thus y may be changed without having to 
redo many of the calculations within the program. 

Intuitively, the idea is to increment x through a range of values, while de- 
termining y at each x according to the criterion mentioned below, until it is 
determined that a solution has been crossed. This determination employs the 
use of a test function which changes sign when a root is crossed in the same 
manner that a single equation in one unknown changes sign as it goes through 
a zero. The task of developing a defining criterion for y and a test function 
for x would be easy i*, for example, v were a strong function of y and u were 
a strong function of x. Then, for each x, y could be chosen such that 
v(x,y) = 0 and, as x is incremented, a solution would be cro c ’ when u[x,y(x)] 
changes sign. Unfortunately, neither u nor v behaves this wa j 

To develop functions that do behave this way, the following procedure was 
developed. First, for each x choose y such that the absolute value of F is 
minimized. This can be done in several ways. The program uses a method that 
always guarantees finding a value if one exists. Essentially, the absolute 
value of F squared and its gradients are computed. The value of y is altered 
in the direction indicated by the gradient until either the gradient changes 
sign or is so close to zero that convergence has been reached. Once the 
gradient changes sign, Muller's method (Ref .28) is used to converge on the 
root. This is essentially a bisection method followed by inverse parabolic 
interpolation. For this searching process, it is not necessary to make the 
convergence criteria very tight, since only rough estimates are eventually 
needed in order to start the two-dimensional secant method described earlier. 
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Now that a criterion for y has been established, it is only necessary to find 
a test function that will chanqe sign when a solution is crossed while incre- 
menting x. Such a function is given by 


uu x + vv x (136) 

This function acts as a very good test function because it represents the co- 
ordinate direction in the u,v plane along which the vector (u,v) changes most 
with x. When this coordinate changes sign as one goes from, say, x-j to 
with y.| and y^ chosen so that the length of the vector (u,v) is minimized, 

then it is very likely that a solution has been crossed. Exceptions to this 

2 2 

rule occur when one is near relative minima of the surrace u + v that are 

not zero. To see this, consider the actual equations that are being solved. 

In order that the vector (u,v) is minimum for each y, it is necessary that 
? 2 

3(u + v )/?y = 0. That is, uu y + vv y = 0. Combining this with the above 
equation, we see we are finding an x and y such that the matrix equation 



(137) 


is satisfied. The matrix is just the transpose of Jacobian of u and v with 
respect to x and y. 


This equation can be satisfied if either u and v are zero, or the Jacobian 
is singular. The Jacobian is necessarily singular at all relative minima of 
the surface u c + v except those at u = v = 0. In order to differentiate 
between those solutions to(137)that are due to singularities of the Jacobian 
ar, those that are due to u and v vanishing, we employ two different tests. 
First of all, when a singularity point is crossed, the determinant of th_- 
Jacobian should charge sign. If this occurs, then the program rejects this 
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point as a possible solution. Sometimes, however, the determinant does not 
change sign because either the convergence criterion used in the searching 
algorithm is too loose or because the singularity has a double root. In 
either case, the procedure is to test the condition number* of the transpose 
of the Jacobian matrix in the region near the suspected solution. If the 
condition number does not exceed a given input limit (e.g. , around 80), 
then the point in question is usually a solution. 

Once it is determined that a potential solution has been crossed between x 1 and x,, 
for example, the procedure is to (a) determine x-, based on the method of 
false position using the test function given in(136),(b) determine y 3 to 
minimize |F|, and (c) use (xj.yj), (x 2 ,y 2 ), and (x 3 ,y 3 ) as the required first 
th.'ee guesses for the two-dimensional secant method. Operationally steps (a) 
and (b) are repeated (at least twice) until the total error, |F|, is less 
than a given fraction of Ja n ^ . This has been found to be necessary since 
occasionally one has to be quite close to a singularity before it can be 
discovered. Repeating steps (a) and (b) will result in a convergence to that 
singularity which will then be discovered by either the determinant changing 
sign or the condition number getting large. However, once the total error 
is reasonably small, the determinant remains of one sign, and the condition 
number is not large, it is almost certain that an actual solution exists 
and therefore the program proceeds to step (c). 

The above procedure has been found to be most satisfactory for the conditions 
tested in this program. The search algorithm described above has several 
salient features. First, as mentioned earlier, the search method takes ad- 
vantage of the fact that "any of the computations in the program are not a 
function of the imaginary part of w, namely y. Specifically, the hydro- 


* The condition number of a matrix. A, is a measure ot how sensitive a 
solution to the system Ab=c is to perturbations in c. It is equal to the 
square root of the ratio of the absolute value of the largest eigenvalue of A’A 
to the smallest eigenvalue of A'A. For singular matrices, the condition 
number is infinite. For matrices that are nearly singular, the condition 
number will be quite large. 
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dynamics of the feed system, the combustion dynamic coupling coefficients, 
and the major terms needed to compute the downstream nozzle admittance are 
all only a function of the real part of the frequency. This allows the 
minimization of |F| with respect to y to proceed with high efficiency. 
Secondly, and more importantly, the procedure has been automated to the 
extent that the user only has to specify a frequency range and a maximum 
number of roots he desires in that range. The algorithm will start at the 
lower end of the frequency range and will increment through it until either 
the maximum number of roots are found or the upper end of the frequency 
range is reached. This is a very powerful property since it does not require 
the user to have a clear knowledge of the location of any of the roots in 
the 10 plane. 
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SECTION IV 
MODEL VERIFICATION 


INTRODUCTION 

The primary objective of Phase III of the OME/Feed System Coupled 
Stability Investigation was to verify the operation and capabilities 
of the Engineering Model, and to compare the model's predictions 
with experimental data from an OMS-type engine with known feed system./ 
engine chugging or buzzing history. In the following sections, the 
experimental hot-fire test data which was used for comparing the model's 
operation and performance are described. Next, the methods of evaluating 
the combustion parameters and Klystron constant appearing in the com- 
bustion dynamics model are presented. Finally, e summary of the Feed 
System Coupled ability Model results from each of the seven model 
verification tests is given. The computed output is compared to ^ihe 
experimental predictions and analyzed in terms of the applicability of 
the available experimental data to the OME system being modeled. 

DESCRIPTION OF EXPERIMENTAL HOT-FIRE TEST DATA 
Background 

Rocketdyne has recently completed an extensive and conorehensi ve program 
for NASA/Johnson Spacecraft Center (NAS9-12802) to determine the feasi- 
bility of and evaluate potential reusable thrust chamber concepts for the 
Space Shuttle Orbit Maneuvering Engine (OME). The program determined the 
applicability of various thrust chamber concepts to the Space Shuttle OME 
application. Feasibility was analytically predicted and experimentally 
demonstrated for the most promising reusable thrust chamber concept. This 
information will support NASA/JSC and shuttle vehicle contractor Orbit 
Maneuvering System (OMS) studies, and provide technical foundation for the 
final definition of the OMS. The technical effort started in June 1972 and 
was completed in February 1975. 

Under Task XII of the SS/OME Reusable Thrust Chamber program a series of 
experimental tests was conducted at White Sands Test Facility to investigate 

^RECEDING PAGE BLANK NOT FILMED 


R-9807/IV-1 



(1) the start, shutdown, and restart characteristics of the integrated 
thrust chamber; and (2) OME thrust chamber operating characteristics 
at very low chamber pressures typical of propellant tank blowdown opera- 
tion and without supplementary boundary layer coolant. Test Sequence 7 
of this task involved a series of tests in which the level of chamber 
pressure was progressively reduced until chugging occurred and included 
the effect of mixture ratio on the minimum pressure level. Discrete 
5-second tests were conducted (as opposed to continuous blowdown tests) 
to minimize facility hyperflow time. Cold fuel (45 F) was used to enhance 
the regenerative coolant safety factor at low chamber presrures. Tests 
7-7 through 7-14 of this series were conducted with unsaturated propellant 
using the facility feed-system configuration which simulates the OMS. 

Tests at 80 ar.d 75 psia nominal chamber pressure resulted in chugging at 
start which damped out durinq the tests. At 65 psia nominal chamber pres- 
sure, the chug persisted throughout the test. Chugqinq was observed to 
occur at frequencie c between 115 and 350 Hz. 

After extensive review of Rocketdyne's hot-fire test data rbtained at 
White Sands Test Facility for the SS/OME Reusable Thrust Chamber Proqram 
and following discussions with the NASA/JSC contract monitor, it was con- 
cluded that the data reported in Task XII (OME Integrated Thrust Chamber 
Tests) would be sufficient and appropriate for use in Phase III model veri- 
fication of the current program. An "Engineering Model Verification Plan" 
(Ref. 29) was subsequently written which summarized all pertinent White 
Sands Test Facility hot-fire test data to be used during model verification 
The report included (1) details of Rocketdyne's thrust chamber and injector 
(like-doublet No. 1), (2) schematics of the fuel and oxidizer feedline con- 
figurations, and (3) steady-state operating data for the seven planned veri 
fication analyses. A summary of the pertinent test hardware, test facility 
description, and hot-fire test data contained in this report is presented 
below. The reader is referred to Ref. 30 for a complete description of all 
test hardware, instrumentation, and experimental results pertaining to 
NAS9- 12802, Task XII. 
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Test Hardware 

The hardware used for the White Sands test program (Task XII) consisted 
of a regeneratively-cooled thrust chamber, full size and truncated radia- 
tion-cooled nozzles, and a like-doublet injector. The injector and 
chamber were designed to closely simulate the thermal and dynamic charac- 
teristics of OME flight-type hardware. All components were bolted together 
and sealed with either metallic or elastomeric 0-rings, as appropriate. 

Table 2 provides a summary of the regeneratively cooled chamber design 
characteristics. The combustion chamber had a length of 0.3/34 m (14.7 in.) 
and a contraction ratio of 2:1 with a throat diameter of 0.1478 m (5.820 in.). 
The expansion area ratio of the regeneratively cooled nozzle was 7:1. A radia- 
tion-cooled Columbian nozzle with an expansion ratio of 9:1 was used for the 
tests at low chamber pressure (Test Sequence 7) to eliminate the chance of 
chamber damage due to high side loads. The inner wall and the lands of the 
chamber were 321 CP.ES, and the channels were closed out with electroformed 
nickel. The thrust chamber was designed for the heat flux profile shown in 
Fig. 13. Channel sizes were such that the minimum safety factor was approxi- 
mately 1.5 at a fuel inlet temperature of 38.1° C (100° F), chamber pressure 
of 82.73 x 10^ N/m^ (120 psia), and propellant mixture ratio of 1.85. The 
coolant jacket itself was flightweight with nickel closeout thicknesses as 
thin as 0.064 m (0.025 in.) at the throat. The fuel inlet manifold was a 
heavyweight configuration to reduce cost, but simulated flight manifold 
volume. The coolant outlet manifold was more critical thermally and repre- 
sented a typical flight design. 

The completed regeneratively-cooled thrust chamber is shown in Fiq. 14. 

Complete details of all instrumentation used durinq Task XII is given in 
Ref. 30. It should be noted that chamber pressure instrumentation was 
not configured to provide high response transient data. However, close 
coupled high-pressure, high-response transducers were used to record the 
transients in fuel and oxidizer injection pressures and fuel coolant 
jacket inlet pressure. Other principal instrumentation involved measurement 
of propellant flows, thrust, chamber skin temperatures, and P c in acoustic 
cavities . 
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TABLE 2 . DEMONSTRATOR THRUST CHAMBER DESIGN CHARACTERISTICS 


COMBUSTOR 

Contraction Ratio 
Length, m (in.) 

Contour 

NOZZLE 

Regen Section Expansion Ratio 
Nozzle Extension Expansion Ratio 
Contour 


2:1 

0.3734 (14.7) 

Tapered from 0.1778 m 
(7 in.) upstream of 
throat 

to 7:1 

7:1 to 72:1 

Flight parabolic 


COOLANT 

Circuit 

Number of Regen Coolant Channels 
Coolant Pressure Drop, N/m^ (ps id) 
Coolant Bulk Temperature, Rise, °C (°F) 
Auxiliary Film Coolant 
Channel Dimensions at throat 

Width, m (in.) 0.0029 (0.114) 

Height.m (in.) 0.0017 (0.068) 

Channel Dimensions near in.iectcr 

Width, m (in.) 0.0029 (0.114) 

Length, m (in.) 0.0011 (0.042) 


Counterflow 

120 

10.34 x 104 (15) 

81.4 (178) 

2.7% Total Propellant 


MATERIALS 

Hot Wall (0.0008 m/0.030 in.) and Lands 
Cold Wall (0.0008 m/0.030 in.) 

Nozzle Extension 


CRES 321 

Electroformed Nickel 
CRES 
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The Injector used was a like-doublet (L/D No. 1), which had 186 elements 
arranged In nine rows. Oxidizer orifice diameters ranged from 0.0008 to 
0.001 m (0.032 to 0.038 in.), while fuel orifice diameters ranged from 
0.0007 to 0.0008 m (0.028 to 0.033 in.). The Injector Included 68 orifices 
(0.0005-m, 0.020-in. diameter) to provide boundary layer coolant amounting 
to 2.7 percent of the total propellant flow at nominal mixture ratio. Ip- 
jector characteristics are summarized in Table 3. 

The injection pairs of like-double No. 1 were designed with a fan impinge- 
ment angle of 60°, an oxidi zer- to-fuei fan inclination angle of 22.5°, and 
an element spacing 0.0073 m (0.289 in.). These parameters are shown 
in Fig. 15. The nominal impingement height of the like-doublet pairs was 
0.0048 m (0.188 in.). Detailed iniector orifice data are presented in 
Table 4. 

Test Facility 

The thrust chamber assembly was tested at the White Sands Test Facility at 
Las Cruces, New Mexico. A simplified feed system schematic of the in- 
stallation is shown in Fig. 16. Fuel (MMH) and oxidizer (NT0) were stored 

3 

and conditioned in 7.57 m (2000 gallon) propellant tanks external to the 
vacuum cell. The propellant was pumped from the external tanks to the two 

3 

0.227 m (60 gallon) tanks inside the vacuum cell simulating the 0MS tank- 
age exits. Line sizes and lengths from the propellant tanks to the 0ME 
interface were configured so as to simulate 0MS ducting. Dimensions of 
the actual fuel and oxidizer feedline configuration used during integrated 
hardware testing are given in Fig. 17 and 18. A flowmeter was located in 
each propellant feedline between the tanks and the engine interface. A 
common pressure source was used to pressurize both the internal anu extern; 
propellant tanks. 

The fuel sides of the two LM descent engine valves were used as the engine 
propellant control valves. Fuel valves were used for both fuel and oxidizer 
sides because these valves contained the actuators and the position indica- 
tors. Each valve was series and parallel redundant including upstream iso- 
lation valves and downstream shutoff valves. ” >itions were measured on one 
of the isolation valves and one of the shutr'f valves for each propellant. 
The valves were located so as to provide a slight positive drain into the 
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TABLE 3. INJECTOR .,/D No. 1 CHARACTERISTICS 


Diameter, m (In.) 

0.2083 '8.200) 

Number of Elements 

186 

Number of Rows 

9 

Type of Elements 

Like Doublet 

Oxidizer Eleme -neter, m (in.) 

0.0008/ / 0.032/^ 

(mini mu>.i/ max i mum) 

0.0010 \ 0.038 / 

Fuel Element Diameter, m ( i .. . ) 

0.0007/ / 0.028/ \ 

(mini r«em/max i mum) 

0.0008 \0.033 / 

Pressure Drof § Nominal Flows 


Oxidizer, N/lli 2 (psi ) 

38.61 x 10* (56) 

Fuel, N/m 2 ipsi) 

42.75 x 10 4 (62) 

Number of Acoustic Cavities* 

8/4 

Mode Suppression 

1 r.i & 3rd 
Tangential , 
1st Radial 


* Cavities formed by chamber and injector 
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Figure 15. Like Doublet No. I Element Configuration 
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TABLE 4. OME LIKE DOUBLE T INJECTOR PARAMETERS, LIKE DOUBLET NO. 1 






INJECTION 

ROW 




ACOUSTIC 
C VITY 
DAM BLC 


1 

2 

3 

4 

5 

6 

7 

8 

9 


.0345 





.0345 

•0321 

.0332 

.0379 


UXIO* UT1 ; 3 Uia 





m 

Fuel Orifice Die 

.0300 

^ 




.0300 

.0279 

l 

e 1*89 

.0330 

.0200 



f" 


No. of Elements 

6 

10 

12 

16 

20 

24 

30 

34 

34 

42 

Element Spacing 

.754 

.688 

.770 

.725 

.697 

.679 

.622 

.618 

.669 

- 

Row Radius 

.720 

1.095 

1.470 

1.845 

2.220 

2.595 

2.970 

3.345 

3.620 

- 

Face Area Fract. 

.0515 

•0513 

.0689 

.0865 

• 1041 

.1216 

.1392 


.2112 

m 

Flow Coverage 

.0321 

• 0534 

.0641 

.0855 

.1069 

.1283 

.1392 


.2207 

m 

Flow/Area 

.629 

1.041 

.931 

.989 

1.027 

1.054 

1.000 


1.045 

m 

BLC % Fuel 










1.3 

Design Delta P 

AA A 





l 



44.8 

Q 


UAitt* aC * 

Fuel 

W • O 

52.8 

1 1 




e 








52# o 

o 

1/2 Thrust/ 
Element 

16.05 

16.02 

16.03 

16.03 

16.04 

! 

16.04 

13.92 

14.98 

19.48 

- 
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Figure 16. Simplified NASA/WSTF Propellant Feed System Schematic 
With Ducting to Simulate QMS 




TANK 
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1 .803 m 
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FUEL 
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Figure 17. Rocketdyne Integrated Hardware Tests, Fue< Feedline Configuration 
Simulated APS POD Lines at White Sands Test Facility (3-74) 
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215° 


Figure ' 8 . Rocketdyne Integrated Hardware Tests, Oxidizer Feedline Configuration, 
Simulated APS POO Lines at White Sands Test Facility (3-74) 





engine Inlets In an attempt to simulate the depletion which would occur 
after shutdown under zero ' g ' conditions. The ducting between the valves 
and the engines was configured to simulate typical line volumes and sizes 
for the flight OME as shown in Fig. 19 . 

Details of the test facility instrumentation and test stand operation are 
given In Ref. 30. 

Hot-Fire Test Data 

As stated earlier, one purpose of Test Sequence 7 was to Investigate OME 
operation at very low chamber pressures typical of propellant tank blowdown 
operation. To prevent excessive side loads which might be associated with 
very low pressure operation, a 9:1 nozzle was used. In addition, the fuel 
temp ature was reduced to approximately 7. 5*C to provide a higher safety 
factor. 

Tes's 7-2 and 7-4 were moderately low-pressure tests conducted with the 
external facility propellant tanks. Tests 7-7 through 7-14 were conducted 
with unsaturated propellant using the facility feed-system configuration 
which simulates the OMS. Based on the stiffness factor, YaP/w (where &P 
is the difference between propellant tank and chamber pressures, and w Is 
the propellant flowrate) the feed system using external oxidizer tanks is 
3% stiffer than the OMS simulated system and the feed system using external 
fuel tanks is 8t stiffer than the corresponding OMS system. 

The range of chamber pressures, propellant flowrates, and mixture ratios 
employed in Tests 7-7 through 7-14 is shown in Table 5 . Measured and 
predicted injector pressure drops are tabulated. The measured Dressure 
drops are the values recorded when no significant oscillations were oc- 
curring during the test. Predicted values are all referred to the measured 
values of Test 7-7 and ratioed according to the sauare of the flowrates. 
di f f a rence between the measured and predicted pressure drops partlcu- 
ar v 1 ower pressures may be the result of subtracting two high pres- 
* . aments (injection and chamber pressures) to obtain a aP although 
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Figure 19. Propellant Inlet Ducting and Valves 




TABLE 5. SS/OME REUSABLE THRUST CHAMBER PROGRAM (NAS9-12802) 

TASK XII - OATA OUMP 
INTEGRATED CHAMBER LOW-PRESSURE TESTS 



7-7 

7-8 

7-9 

7-10 

7- IDA 

7-11 

7-12 

7-14 

ChAMBLR PR. N/M 2 X I0 4 (PSlA) 

80.37 

57.92 

52.40 

46.88 

50.33 

46.50 

46.33 

36.88 


d? 6 ) 

(04) 

(76) 

( 68 ) 

(73) 

( 66 ) 

( 6 £ v 

(126) 

HIXTURI RATIO, U/l 
( LObRATL . KG/SLC (LU/SLL) 

1 .68 

1.7 

1 . 8 . 

1.8 

1.6 

1 .6 

1 

1 .66 

o\w 

5.62 

3.76 

3 . 3*9 

2.33 

3.08 

2.73 

2 . <58 

6.53 


(12.4) 

(0.29) 

(7.41) 

(6.34) 

(6.73) 

( 6 . 02 ) 

(6 .67) 

( 12 . 2 ) 

FULL 

3.38 

2.14 

1.36 

1.58 

2.06 

1.7? 

1.54 

3.36 

MLASURLU inj p.* 

:./m 2 \ lo 4 (PSD oxiu 

(7.39) 

(4.72) 

(4.10) 

(3.43) 

(4.53) 

(3.79) 

(3.40) 

(7.39) 

38.61 

13.10 

7.58 





36.lt) 

(56) 

(19) 

on 

-- 




(51) 

ULL 

39.30 

.4 I 

6.39 





35.85 


(^7) 

x 18) 

(1C) 



-- 

-- 

:s 2 ) 

PRLDICTLO Iw. P , 









;./m 2 X 13 J (PL I) OX 10 

38. tl 

16.66 

13.79 

in. 34 

11.72 

3.96 

P .03 

37.23 

(56**) 

( 24) 

( 20 ) 

(16) 

(17) 

(13) 

U 6 ) 

(54) 

1 UlL 

39.30 1 

15.36 

11.72 

3.27 

14.43 

10.34 

8.27 

39.30 


(57** ) 

(23) 

(17) 

(ID 

( 21 ) 

(lb) 

( 12 ) 

! j7 ) 

PRLDICTLO P/ Pc OX I lx 

9.44 

3.2* 

0.26 

0.22 

o. 2 J 

0.20 

0 24 

0.42 

FuLL 

0.46 

0.23 

0.22 

n 19 

i n .29 

0.23 

C.1S 

n 43 

FPLIJULUCY. Hi 
PK/PK AtIPLITUOL, 

STABLE 

STABLE 

STABLE 

246 

360 

280 

255 

STA3LE 

:,/m 2 x io 4 (ps I) ox id 


4.82- 

6.39 

49.64 

6.21 - 

23.44 

33.no 

S.27 


11.72 



! 16.54 






(7-17) 

( 10 ) 

(72) 

;°-24) 

(34) 

(48) 

( 12 ; 

FotL 


4.32 

4.82 

24.82 

6.89 

16.65 

20.27 

6.39 


(7) 

(7) 

(36) 

( 10 ) 

i 

(24) 

(41) 

( 10 ; 


•DURING STACIE PA '7 Or TEST 
••KLFLKUU 
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the discrepancy is greater than would be expected from instrumentation 
inaccuracies. Values of injector AP divided by P £ are also tabulated 
based on the predicted pressure drops. Tests were conducted with values 
as low as approximately 0.2 on each side of the injector. 

Chugging was not evident as the chamber pressure was reduced from test 
to test until Test 7-8 where a chamber pressure of 57.92 x 10^ N/m^ 

(84 psia) was obtained for steady-state. On this test, chugging occurred 
for approximately 0.6 seconds after ignition at a frequency of about 300 
Hz. As the chamber pressure was reduced on the successive tests, the 
duration of the chugging increased until at a chamber pressure of 46.88 x 
10^ N/m^ (68 psia) chugging continued throughout the test. When the cham- 
ber pressure was increased on the next test to 50.33 x 10^ N/m^ (73 psia) 
and the mixture ratio was decreased to 1.5, <ontinuous chugging occurred 
for 1.6 seconds but continued sporadically throughout the remaining dura- 
tion of the test. Variation of the mixture ratio to 1.6 and 1.9 at 
46.19 x 10^ N/m^ ( 67 psia) chamber pressure on the next two tests resulted 
in small changes in the chug amplitude. Chug frequencies of 200-300 Hz 
were noted. Oscillations at two frequencies were indicated during the 
early portions of some of the tests; however, the lower frequency persisted 
for less than 1 second after start. 

Complete summaries of all steady-state performance data (average measured 
values and calculated values) compiled for Test Sequence 7, Runs 8-14 are 
presented in Ref. 29. 

Application of Experimental Data to Model Verification 
The White Sands test data presented in the previous section are sufficient 
to provide initial verification of the engineering model in a manner which 
is consistent with the modtl' c anticipated normal usaqe. A total of seven 
verification analyses were evaluated based upon the integrated chamber low 
pressure tests s’.vnarized previously in Table 5 . Tests 7-8, 7-9, and 
7-14 were used to represent examples of staole system benavior while tests 
7-10, 7-11, and 7-12 provided examples of unstable engine operating condi- 
tions. The low amplitude chugqlng behavior encountered In test 7-1 CA was 
planned to serve as an illustration of marginally stable operating behavior. 
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MODEL INPUT DATA 

Each of the seven verification analyses reoulred the same basic Input 
data to the computer program, namely: 

1. Feed system configuration (Figs. 17 and 18 ). 

2. Injector characteristics (Fig . 15, Tables 3 and 4). 

3. Chamber geometry (Fig. 14 , Table 2 ) . 

4. Steady-state engine operating conditions ( Re f- 29 Appendix). 
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In order to determine the steady-state distributed enerqy release functions 
required by the stability model, an existing Distributed Energy Release 
combustion model (JANNAF computer program DER) was executed for each veri- 
fication case. In all cases, the DER calculated performance was less than 
the measured engine performance. This difference between the calculated and 
experimental performance should not affect the stability calculations since 
performance is dependent upon the largest droplet diameter groups while 
stability is mainly dependent upon the smaller droplet groups. 


Output from the DER computer calculations was curve- fitted with the steady- 
state equations employed by the stability model in order to obtain important 
combustion-related parameters. One of these parameters Is the vaporization 
"time delay''. Integrating the steady-state vaporization expressions qlven by 
Eq. 64 from the start plane for vaporization (X Q ) to any given vapori- 
zation plane (X) yields: 


In (% unburned propellant) 


= ln(100) 



(138) 


A plot of this equation along with DER calculated values for the fuel and 
o v idizer for test IHT 1-7-14 is shown in Fig. 20. A straight line was drawn 
through eacii curve for X < 0.0635 m (2.5 in.), since the vaporization rates 
corresponding to these distances are much greater than those for X > 0.0635 m 
(2.5 in.). Obtaining t V from the slope of each curve and using the fact that 
V equals the average injection velocity enables the vaporization "time delays" 
to be calculated. Methods for obtaining other input data required by the 
stability model from the DER program output are outlined in Ref. 3. 


Initial model verification calculations indicated that the following terms 
must be neglected In order to obtain agreement between calculated and ex- 
perimental results: 
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1. Changes In the vaporization time delay due to changes In 

mixture ratio must be neglected, l.e., = 0. 

2. Changes in the nozzle admittance due to mixture ratio chanqes 

3c* 

must be neglected, l.e., “ 0. 

3. The oscillatory relative velocity between the gas and the droplets 
must be neglected, i.e., x DRAfi * 0. 

Since the vaporization "time delay" is actually temperature dependent instead 
of mixture ratio dependent, even though the temperature .s a function of 
the mixture ratio, the stability model should be modified to reflect this 
dependence. Model calculations at the end of the program indicated that if 
this change from mixture ratio to temperature dependence is made in the model, 
agreement between calculated and experimental results can be maintained while 
ir,. luring (|y) and effects. Problems encountered when the oscillatory 

relative velocity effects are included in the model calculations have not been 
isolated but continuing effort will be directed towards including this im- 
portant coupling term. 

DETERMINATION OF KLYSTRON DISTANCE 

One of the important input parameters required by the stability model is 
the Klystron distance. Since experimentally determined /alues for the 
Klystron distance for the OME system configuration are not presently avail- 
able, and in lieu of actual combustion characterization test data, tte Klystron 
distances were determined based on IHT 1-7-10 and 1-7-12 test data. The 
fuel and oxidizer Klystron distances were assumed to be equal and independent 
of engine operating conditions. These assumptions are probably invalid but 
were made for lack of actual combustion characterization test data. 

Shown in Fig. 21 is the decrement for tests 10 and 12 as a function of the 
Klystron distance. The decrement (x) is defined as minus the imaginary part 
of the complex frequency divided by the real part of the complex frequency 
and is a measure of how fast the wave will grow or decay. A positive value 
for the decrement indicates a damping or stable behavior, whereas a negative 
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Figure 21. 
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Determination of Klystron Distance Based upon W5TF Tests 
IHT 1-7-12 AND IHT 1-7-10 
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value indicates an increasing amplitude or unstable behavior. Marginal 
stability or instability corresponds to a decrement value near zero. 

Since tests 10 and 12 were experimentally unstable. Fig. 21 yields the 
result that the Klystron distance must be qreater than 0.0495 m (1.95 in.). 
For model verification, a Klystron distance of 0.0508 m (2 in.) was selected 
based on the information presented in Fig. 21. 

COMPARISON BETWEEN CALCULATED AND 
EXPERIMENTAL RESULTS 

Execution of the complete engineering model results in the prediction of 
the oscillatory frequency and information related to the stability or damp- 
ing coefficient (decrement) of that frequency for the pressure-fed propul- 
sion system. For each given test case, the frequency and stability predicted 
by the engineering model can be compared with the experimental data. Shown 
in Table 6 are the tabulated experimental data and calculated results. 
Comparison between calculated and experimental results indicate that the 
model is in agreement with the experimental data in predicting stable or 
unstable operation for all tests except IHT 1-7-11. Also, the value for 
the decrement for test IH~ 1-7-10A should be closer to the neutral stability 
value of zero. Comparison of the experimental and calculated frequencies 
for the unstable tests yields the result that the calculated Instability 
frequencies are much lower than the experimental values. This discre- 
pancy in frequency could be due to assumptions in modeling the feed system 
or to omitting certain combustion effects. 

Figure 22 presents the experimental and calculated results as functions 
of the fuel and oxidizer injector pressure drops. The experimental data 
can be separated into stable and unstable regions, i.e., if the fuel and/or 
oxidizer delta pressure is below a certain value, the engine was experi- 
mentally unstable. Based on this fact and the calculated results presented 
in TaDle 6 , Fig. 22 indicates that the model calculations are not sensitive 
enough to the oxidizer Injector pressure loss. Based on tests 9 and 10 
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TABLE 6 . COMPARISON OF CALCULATED 
AND EXPERIMENTAL RESULTS 



EXPERIMENTAL 
FREQUENCY (Hz) 

CALCULATED 

TEST NO. 

FREQUENCY 

DECREMENT 

8 

Stabl e 

216.4 

0.1370 

9 

Stable 

209.6 

0.0417 

10 

245 

207.1 

-0.0078 



281.0 

0.1024 

10A 

Marqinallv 

152.1 

0.1431 


Stable (350) 

175.8 



211.3 




247.3 

0.1535 

11 

280 

198.6 

0.1748 



219.8 

0.1598 



261.4 

0.1159 



293.2 

0.0695 

12 

255 

174.3 

0.1901 



212.8 

-0.0533 



280.4 

0.09566 

14 

Stable 

264.8 

0.2613 



314.8 

0.2020 
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Figure 22. Comparison of Feed System Coupled Stability Model 
Calculations and WSTF Experimental Data 




(or 12), the model predicts the correct trend with respect to fuel delta 
pressure but examination of tests 10A and 11 shows that the model does not 
predict the correct trend with respect to oxidizer delta pressure. This 
insensitivity to oxidizer delta pressure could be due to the particular 
hydrodynamic modeling of the oxidizer system or due to using the same 
Klystron distance for the oxidizer and fuel. It should be noted that the 
Klystron distance was determined based on tests 10 and 12 which are con- 
trol lec by the fuel system. Also, the omission of certain combustion effects 
could, and probably are, affecting the model predictions. Final determination 
of the models inability to predict correctly the experimental results of tests 
11 and 10A cannot be made until after confcustion characterization studies 
have been completed. A limited number of calculations have shown that the 
model is quite sensitive to changes in the hydrodynamic system and certain 
combustion parameters. 
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SECTION V 


CONCLUSIONS AND RECOMMENDATIONS 

The Feed System Coupled Stability Model (FSCSM) has been formulated* 
checked out and successfully executed. It has been shown to be 
capable of analyzing feed system/engine coupled instabilities in 
pressure-fed, storable propellant, propulsion systems over a 
frequency range of 10 to 1000 Hz. The model has been constructed in 
a generalized manner In order that it may be readily applicable to 
present and future engine and propulsion systems. 

3ased upon the limited number of model verification analyses 
conducted thus far using the OME system as a working basis, the 
model appears to predict correctly the stability or instability of 
a given engine system. As expected by the theoretical analysis, 
the Klystron constant was determined to have a significant effect 
on the stability prediction. With respect to the predicted unstable 
tests, the frequencies calculated by the model were founc to be much 
lower than the experimentally determined frequencies (213 Hz vs. 255 Hz 
and 207 Hz vs. 245 Hz). V.ese results are due in part to the assumptions 
used in modeling the feed system as well as to the omission of 
certain combustion effects. The calculated frequencies were found 
to be highly sensitive to the feed system input data, particularly 
to the lengths of the various line segments and to the injector 
manifold volume. In addition, the value of the decrement was directly 
related to the magnitude of the Injector HP. Comparison of the 
available experimental and calculated results Indicated that the 
model was not sensitive enough to the oxidizer Injector pressure loss. 

A number of areas exist where modifications or possible Improvements 
to FSCSM can be made. These are listed below: 
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1. Input data required by the model should be simplified. This 
could be accomplished In two steps. First, modify the output 
of the steady-state DER program to calculate and print out In 
a logical sequence the pertinent data to be used as Input to 
the combustion and chamber dynamics subprograms. Next, change 
the hydrodynamics subprogram to accept line segment AP's as 
Input data and calculate values of resistances and acoustic 
velocities Internally. 

2. More accurate feed system formulations should be Investigated 
to Include the effects of various nonl Inearl ties. This would 
Involve the use of the complete continuity and momentum equations 
Instead of the simplified "waterhammer" description. 

3. Improvements in the combustion and chamber dynamics model 
formulations should be studied. First, the velocity coupling 
terms should be reviewed. These were neglected during model 
verification due to their extreme sensitivity in affecting the 
model results. In addition, more realistic methods of evaluating 
the effect of mixture ratio on nozzle admittance should be 
determined. Improved methods of calculating the vaporization 
time delays as functions of pressure, temperature, mixture 
ratio, etc., should also be Investigated. 

Inclusion of the full set of spray equations (continuity, 
momentum, etc.) would enable the Klystron constant to be 
determined more readily as a function of the propellant properties 
and/or Injector geometry. This would also yield more realistic 
expressions from >ich the vaporization rates and droplet 
velocities could be calculated. Although this would tend to 
Increase the complexity of the chamber dynamics equations, the 
solution method for determining the complex frequency and 
decrement (stability rating) would not be affected. 


R-9807/V-2 



Finally, to Incorporate certain types of non- linearities into 
the combustion dynamics formulation, the combustion coefficients 
should be made functions of the local peak-to-peak pressure 
amplitude. Then, Investigation of stability or Instability for 
a finite disturbance could be made and the actual peak-to-peak 
amplitudes determined for unstable disturbances. 

4. Actual OME subscale test data should be obtained yielding 
combustion characterization test parameters and detailed feed 
system Information to provide a more rigorous verification of 
the FSCSM. The experimental test program should be structured 
Into two phases. The first phase would be used to check out 

the hydrodynamics model formulation. It would Involve perturbing 
the pressure In the test chamber under non-combustion conditions 
and measuring local oscillatory pressures and flowrates at 
selected points along the length of the teed system. The 
second phase would ? *'• verify the combustion and chamber 

dynamics models. c,.*M'.1st of pulsing the feed system 

at a known frequer..^ <nd measuring the combustor response. 

5. Additional areas exist where improvements in the model could 
aid In better understanding of typical feed sy'tem coupled 
instabilities. These Include: 

a. Multi-dimensional chamber and manifold descriptions 

b. Couplir.g of pressure response with Injector passage 
resonance 

A multi -dimensional model would be very useful in analyzing instabilities 
where response In the chamber is not uniform across the Injector face. 
This would Include analyzing the effect of feed system coupling on 
chamber radial and tangential modes of Instability as well as 
Including the proper feed system effect In the design and analysis of 
acoustic cavities. Previously, feed system contributions to chamber 
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resonant modes of Instability have usually been Ignored. However, 
some flow oscillations can occur across the Injector face for even 
very high frequency oscillations. In some cases, feed system 
modifications may be the method used to stabilize this type of 
Instability. In other cases accurately Including the feed system 
effect could make design of damping devices more reliable. 

Another type of instability that occasionally occurs Is a feed system 
resonance involving only the injector flow path coupling with some 
pressure response feedback. These instabilities are generally high 
frequency (observed Instabilities have been from 4400 to 14,000 
Hz). The injector resonant response characteristic can be accurately 
modeled, but the feedback loop involving chamber combustion or some 
other generation of pressure oscillations has not been simulated. 

This makes evaluation of various modifications somewhat qualitative. 

If the pressure feedback loop could be accurately simulated, the 
determination of required system modifications would be more reliable. 
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APPENDIX A 


FEED SYSTEM HYDRODYNAMICS MODEL EQUATIONS 

In this appendix, the system of 57 equ. 'Ions describing the generalized feed 
system shown In Fig. 7 is presented. The equations are given In the Itiear- 
Ized, La Place transformed format required by the frequency response sub- 
routine. The derivation of the equations Is based upon the mathematical for- 
mulation outlined In Section II. 
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+ (Rj + A]) Wj ♦ [Pj + (Rj • A]) ft,] e * 0 

(A-l) 
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P 6 + ^**9 + a 9 ^ *12 ” ^ P 5 + a 9 * 11 ^ e = 0 


P 6 " a 10 *12 ~ f P 7 + ^ R 10 " a 10 J * 13 ^ e 


' T 10 s 


T 10 s 
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APPENDIX B 


COMBUSTION MODEL EQUATIONS 

KLYSTRON EFFECT 

The dynamics of the liquid propellant jet from the injector face to any 
location in the chamber are described by the continuity and momentum 
equations: 


at ( A j°j) * h ( A j p j v j)* 0 

(B-l ) 

W (Vj“i) + £ ( A j°j v j) ' - A j& 

(B-2) 

Assuming 

p. = constant 

J 

(B-3) 

*2.= n 

(B-4) 

9x 

$ = <t» * $ U any variable). 

(B-5) 

where 

* f (x) (time average value) 

(B-6) 

$ = $'e~ lut (oscillatory value) 

(B-7) 

*'« g(x). 

(B-8) 

the time average continuity and momentum equations become: 

±(¥j )*° 

(B-9) 

<J* 

M 

O 

(B-10) 

or 

Vj * constant 

(B-l 1 ) 

ITj * constant 

(B-l 2) 
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The first-order oscillatory continuity and momentum equations can be written 
as 


a, 

4 + IT. 


3v, _ 3A. 

— J- + v, 

3x j 3 X 


(B-13) 




0 


Combining Eqs. B-7 and B-U yields 
- dv i 

* "s -af 3 0 

which can be Integrated to yield 



the oscillatory jet velocity can be written as 



Substituting Eqs. 3-7 and B-16 Into Eq. B-13 yields 
Aj(-lw) + (vj)^ e ia *^ v j 

J 

= 0 



(B-14) 


(B-15) 


(B-16) 


(B— 1 7) 


(B-18) 


(B-19) 
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Integrating Eq. B-19 and employing Eq. B-17 yields the jet area 



The mass flowrate of the jet Is 

"j • Vj v j 

Therefore, the oscillatory flowrate Is 


i.i. 




Combining Eqs. B-18 , B-20, and B-22 gives 



(B-20) 


(B-21) 


(B-22) 


(B-23) 


DROPLET VAPORIZATION 

For the fuel or oxidizer spray, the droplet continuity equation can be 
written as 


& < A Vk> * - A N k "»ap k <»-»> 

and the vaporization rate Is (Ref. 12) 

1 " D k kf k N “ H k Z|< 

"vap k i— . (B-25) 

where p k Is the spray density (mass of spray pet unit chamber volume), N^ Is 
the number of droplets per unit chamber volume, and 
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(B-26) 


c Nu pi* D 
p v, "k v k *k 


: — - In /_e_\ 

f k % k « \ l P ' P v t J 


Noting that 


ffD k 

p k = N k m k * N k p t k T - 


(B-27) 


the droplet number flowrate can be written as 

A v. p. 

\ * »k A \ ' -nr 1 


(B-28) 


Therefore, Eq. B-24 can be written as 


Jr ("kV * - £ *w k 


(B— 29) 


For steady-state combustion models, the preceding equation (along with Lq. B-25 ) 
Is numerically Integrated allowing the droplet diameter, D k , to vvy along 
the length of the combustor and maintaining constant droplet number flowrate 
(N k ). Combs' (Ref. 20) has shown that changing from a variable droplet 
diameter to a variable droplet number flowrate yields approximately the same 
results for steady-state vaporization. Therefore, In order to simplify the 
Integration for stability analysis, the droplet diameter was held constant 
and the droplet number flowrate was assumed to vary. 


S d 

k 


y. %p k 


^ \ % 


Summing Eq. B-24 over all fuel or oxidizer droplet size groups yields 

IH, 

k 

which can be written as 


* <Vk> * - A * p k sr ’ -4“T“Z c, 


(B-30) 


K k X 


Ap« 


<nr (Ap s v s } * t 7 “ " A m vap. • 


s vap s 


(B-31 ) 


R-98D7/U-4 



where 


■2 


5 "T P k 


Vk> 


, .-LX 

, , y p k (6)Z k k f k %, 

T s " °S k p. 0? c 


*k k x 


Letting Z. , k- , c n be Independent of k and assuming 
* T k p v 


T s = 


yields: 


Pfl c p D 2 
K v. s 


T s = mr kI ""Ru H “ * 

s s 


where 


D s • ("s’lnj 



From Eq. B-31 


d(Ap $ * s ) 

TvT 


dx 

T s v s 


Integrating Eq. b-38 between x Q and x 


Ap s v s * (m s ) x-x 0 exp 


' / x _dx_ 

■ >x„ Vs 


(B-32) 

(B-33) 

(B-34) 

(B-35) 

(B-36) 

(B-37) 

(B-38) 

(B-39) 
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Substituting Eq. B-39 Into Eq.B-31 yields the fuel or oxidizer spray vapori- 
zation rate: 


KW 

%>, ^vT exp 


t _dx 

r T_V 


OSS 


(B-40) 


Using perturbation techniques, the time average vaporization rate can be 
wrl tten as 


m 


vap. 


( "s>X-X e 

wr 


exp 


(x-x 0 ) 

Vs 


(B-41) 


and the oscillatory vaporization rate can be written as 


m vap s "V s ,. 5 -x=x 0 


/ ( " s, x-k ~s \ 



\ 


+ j i + J-) — 

^ V s] Vs, 


(B-42) 


Assuml ng 


v s ^ ^ v s^x=x. 


(B-43) 


yields 


v_ v « 

JL = (_L ) 
v s % ; x=x k 


(B-44) 


Letting p ff , c , and k- be constant, the oscillatory time delay can there- 
R s p v $ T s 

fore be expressed as 


20s 

0. 


h. 

I. 


Nu, 


Nu, 


H e . , 3t s * MR 
1 +( W } f7 
H. 


(B-45) 
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(B-48) 

(B-49) 

(B-50) 

(B-51) 

(B-52) 

(B-53) 
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Hence, 


R„ v'(-W) •— ^-(1-R„ ) 


T dr«g. 


or 


l+1u) T 


drag. 


1 * ( “ T drag k >' 


(B-54) 


'P-55) 


NUSSELT NUMBER 

It may be observed that one of the dominant termv In the oscill? time de- 
lay Is the Nusselt number. For longitudinal nodes., the Nusselt number is 
(Ref. 21) 


Nu h s 2.0 + 0.6 Pr 


1/3 


PDl 

T- ! v - v rI 


1/2 


(B-56) 


Letting 


1/2 


Fp *(-&-) . F„ 

P 


I v-v. | 


“ca w 


1/2 


where 


AM = 


l!^kl 

c 


-•steady state 
the Nusselt number can be written as 


Nu^ * 2.0 + 0.6Pr 


1/3 


pff 

- 1 CAM 


1 / 


!/2 


F F ( -i 

p M 


(B-57) 


(B-58) 


(B-59) 


Expanding the preceding equation Into time average and 'dilatory parts 
yields P 

PjLk cAM 

L u 


flu H * 2.0 + 0.6 Pr 1/3 


1/2 


r p f v 


>.l 


h f\ ' 

-1/ 

r 1 
+ _* 



i 

) r 

p 

r v J 


(B-60) 


(B-61) 


R-9807/B-8 



From the definition of F and the droplet velocity equations: 

(r n 2 i V /4 

F = { v - v k + v - R v _ ~ Z ~~2 > 
v L K v k -J (v-vj j 
- - 1/2 * ' 


h * (|) 


ar*! 


(B-62) 

(B-63) 


Letting 



the following equation can be determined: 



(B-64) 


(B-65) 


CB- .6) 


For small perturbations In the pressure, the linear response factors a. a the time 
average values are 



(B-67) 


(B-68) 
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OSCILLATORY DROPLET TEMPERATURE 

The oscillatory combustion time delay requires the evaluation of the heat 
transfer blockage term which is related to the combustion gas and liquid vapor 
properties. Since the vapor pressure at the droplet surface Is related to 
the droplet temperature, the blockage term also depends on the oscillatory 
droplet surface temperature. For a single droplet, the temperature inside 
the droplet can be described by the energy equation 

w * 7 & (Aff &]♦ H "- 69) 

where k eff is the effective liquid thermal conductivity. For zero recirculation 
inside the droplet, k^ is the standard liquid thermal conductivity. For 
infinite recirculation. k g ^ is also infinite. 

Since 


pc p T « Pc v T + p 


and assuming 


P = constant, c y * constant, k e ^ * constant. 


eff 

Pc.. 


the energy equation can be written as 


3I = ® _L f r 2 JT_] 

at 2 Jr L 9r J 


The oscillatory part of the preceding equation is 

T ' <-«■ 7 3Fp "'] • 


(B-70) 


(B-71) 

(B-72) 


(B-73) 


(B-74) 
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where 


T . T + r * T ♦ T' e‘ 1ut 


(B-7S) 


and the solution Is 
T* * (T 1 ), 


'r=r. 




/-iw' 
sinh / a 

« ~ 7 ¥ " 


(B-76) 


S‘i 


Differentiating the preceding equation and evaluating at the droplet surface 
yields 


+<£).■ 4 )* “ 


r 

a r s) 


Since the heat transfer rate to the droplet surface is 

%' [ 4 w s k eff(^rJ= \ T i • 

the response factor can be evaluated as: 


(B-77) 


(B-78) 


Rj s = ^ r s keff ^ 


-i ♦ ctnh (/^ r ) 1 

r s a ' ® » 1 J 


HEAT TRANSFER BLOCKAGE TERM 

From the analysis of the gas film surrounding a vaporizing droplet, the droplet 
heating rate can be expressed as (Ref. 12): 


Q s = Zk f Nu h 


<™s> . "a. 
Lle^D C P 


IT D. 


(B-80) 


V f J 


where the heat transfer blockage term is 

A 


z = 




In 


k^Nu f) R T^ 




(B-81) 


R-9807/B-11 



Assuming that the tine average droplet temperature Is at "wet bulb" condi- 
tions [(3T s /3t = 0] and the following variables are constant 


pMW v - v 
c_ , f v f 
\ ' k f T* t f " 


5 

Nu, 


the oscillatory heat transfer blockage term can be written as 



the oscillatory droplet heating rate can be written as 


CB— 82) 


(B-83) 


(B-84) 


(B— 85) 


(B-86) 


q; * 


<"° s > 


2 k f "“h L. 

2 ,> L 


(e z -1) 


AH - 

r + -JS& e * z’ 
S V. 


+ K ( e z . ! ) Is 

C P ' T 

p v f s 

The oscillatory heat transfer blockage term then becomes 

v * a (t — - 
2 T s 1 


(B-87) 


(B— 88) 
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Combining the preceding equations with the heat transfer rate given in the 
preceding section and letting 


the heat transfer blockage term becomes: 


(B 89> 





T 

e 





(B-90) 


2o» k c </-1>T s 
v k 

Zk f Nu h 



Combining the preceding «."■ ess ions, and expressions presented in the model 
formulation section, with the general vaporization rate expression (p. 11-25) 
yields the combustion coefficients: 


c 




+ 




c 
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IF 93' 

(B-94) 

(B-95) 

(B-96) 

(B-97) 

(B-98) 

15 = 0 

,3 S (B-99) 


where the subscript s denotes the fuel or oxidizer and 



(B-100) 



(B-101) 


(B-102) 


It should be noted that the oscillatory heat blockage term has been 
neglected based on work presented in Refs. 22 and 23 which indicates 
that this term is not important at low frequencies. 
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APPENDIX C 


CHAMBER MODEL EQUATIONS 

In this appendix, chamber model equations are derived from basic flow 
equations. Complete derivation of the basic equations Is presented In 
Ref. 31. Assumptions used In the derivation of the basic equations 
are: (1) Ideal gas flow Is valid state equation; (2) dilute sprays 

occupy a negligible fraction of chamber volume; (3) the spray can be 
represented by a finite number of dropslze groups; each dropslze group 
contains a large number of locally Identical drops; and, each size group 
constitutes a separate liquid phase and exchange terms between liquid 
phases are not Included; (4) drag contributes only kinetic energy 
to the spray energy equation; (5) secondary "shear" breakup of drops Is 
not Included; (6) negligible coupling between diffusion and thermal 
gradients; and (7) no body forces. 

The following equations can be formulated for the gas phase: 


Gas Continuity 

-IS > 

n j J 

Gas Momentum 


a 

St 


(P0) + V- (py; 


U) = -Vp + V- T 





Equation of State 


p * pRT 


Shear Stress 

t * -v eff [* + (iff)* - | (v u) fj 


(C-l) 


(C -23 


(C-3) 


(C-4) 
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Gas Energy 


d 

3t 




= - 7. q ♦ 7* (ut) + || 


+ 



n jn 

• V “j 



Gas Mixture Ratio 
1^- (pMR) + 7* (pu MR) 


- P ?L*r T? 2 MR - 2|v MRj 2 1 

eff L nrrr J 


- (VMR) • V(p $ eff ) = 



Heat Transfer Rate 

S *- k eff ,T 'I (cS eff> V»1 
1 


Drag Force 




(C-5) 


(C-6) 


(C-7) 


(C-8) 
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Assuming 


p eff ' °* k eff ~ °* ^eff ' °* 
fJ - 0, 0 ]~ 0, uj = u. 

The preceding equations can be simplified for longitudinal modes to*. 


CC-9) 


Gas Continuity 

A ft * !? (A ’ >v) = A ["vap ox "WfJ 


Gas Momentum 

A ft <‘ >v) *rr (Ai>v2) = - A « + A C\ap ox + *vap fu ^l v 


(C-10) 


(C-11) 


Gas Energy 

A It [ p (h + + h C Ap v (h + irll 

* A H + A L"W 0X (h o* * IT* (c -‘ 2) 

+ "yap fll (h fu*r ) ] 

Equation of State 

p * pRT (C-13) 

Gas Mixture Ratio 

A (p MR) + b (Ap v MR) - A(2 MR + 1) (m yap ) - A(MR) 2 (m^ ) (C-14) 
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Assuming 


h ■ rt * (\ e ,) 0 * (fV 1 , lm ■ %>■ 

(C-1S) 

Wj, ' (v, - 1) * * nd 

(C-16) 

R ■ R * + <!!*>,, < mr -"v. 

(C- 17) 


and substituting these relations Into the preceding longitudinal equations, 
yields the following modified longitudinal equations: 


Gas Continuity 

» H ♦ £ (A > »> ' A <"v, Pox * 


Gas Momentum 


av 


jv + »£. = 


p at + pV ax ax 


Equation of State 


p * » T [V ♦ (y&r> < MR • M V] 



Gas Energy 

A $ + A * Hr + v p n (Av) * 



Gas Mixture Ratio 

P !?♦»»!? ■ < w * ’> [%p ox • ™v,r f J 


(C- 16) 


(C-19) 


(C- 20) 


(C- 21) 


(C- 22) 
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Because of the complexity In solving nonlinear partial differential 
equations perturbation techniques were used to simplify the 
governing dynamic equations. Assuming 

♦ = ♦ + «t> (♦ any variable). 


where 

♦ * f (x) 


(C-24) 


and 

* * g (x, t). 


(C-25) 


each perturbation quantity was taken to be of order (e), where ( e) Is a 
small ordering parameter that Is a measure of the wave amplitude. The 
perturbation expressions for each of the independent variables were 
substituted back into the nonlinear partial differential equations, where 
all terms of ...der (e ) or higher are neglected. The resulting time-averaged 
equations can be solved for the time-averaged variables and the oscillatory 
equations can be solved by assuming solutions of the form 


♦ 


4> ' e 


-io>t 


(C-26) 


where <t>' * f (x) and <*> is the complex frequency. The resulting equations 
form a system of ordinary differential equations in terms of the variables 
4' s and can be numerically Integrated by employing boundary conditions 
and iteration techniques. 


Following this approach the perturbation equations were 

expressed as: 

p * p [_] + p' e' iwt ] 

(C-27) 

v = v” + Cg v' e"* ut 

(C- 28) 

T • T [l ♦ T' e-' 1 *] 

(C - 29) 
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(C-30) 


■ p jj + p' e' 1 "^] 
MR * MR + MR' e" 1ut 


= m i 
OX r ox 


"’vaPnv 'vap rtv vap 


+ m' 


ox 


%» Pfu “ \a Pfu + m vap fu 


■1o)t 


•lut 


The time-averaged equations were determined to be: 


Gas Continuity 

lj_ (A p * A + "U* ) 


vap ox vap fu 


Gas Momentum 

£♦£-» 

Equation of State 

p = p R «f T E 1 + ir 


(MIT - MR, 


■> fi 

Gas Energy 

A * £ + '»' ST lA *> * 


■>] 


< Y o " 11 A j "Vap ox [ 4h ox ' 12 m * 

"vap fu K + *5^^ (Wr)2 l 


+ m. 


Gas Mixture Ratio 


- - d HE" 
p v W 


(fiff + 1) Jj 


m - WIT m J 


vap 


ox 


vap fu- 
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(C-32) 

(C-33) 


(C-34) 


(C-35) 


(C-36) 


(C-37) 


(C-38) 
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The oscillatory equations are given by: 


Gas Continuity 




ox 


vaD^„ vap fu ' 


p c a 


Gas Momentum 


v . (liS) + *1 gv + _». dvl . _pj_ d£ 

' c ' c, ax c. dx - " z ax 

0 0 0 p C 

jy 


P c* 


p C. 


Equation of State 

P 1 * P* * T' ♦ V- ^ *' 

D » 

Gas Energy 

p' ^ i] 

+ ii #+ y f gyl + vL 

— dx O' Ldx X dx 


P 

c * 


] 


+ rsr ai »*>■■? 


(y~- 1) 


P c. 


"v»P 0x [ Sh ox 


(C-39) 


(C-40) 


(C-41) 


(C-42) 
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' 2 MR ’ 1 >] ti v.P fu [ 4h fu 
+ ( im>^ ( *> 2 ] - 2 %ap ox (sm'j, * 

* 2 *m u m < 5 *V MR ' J 

Gas fixture Ratio 



The boundary conditions for the steady-state differential eauatlons are 



"V m x.o (,f Vo ' °> 

* < S U (,f Vo ' °> 


2 

Assuming small Mach numbers, i.e., M « 1, the steady-state differential 
equations can be integrated between the start plane for vaporization (x Q ) 
and any location (x) to yield 
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(C-45) 


p * constant = 



♦ * e 


-(x-x 0 )/t s v s 


_ (Av), 

T 


v = + H£ 


(Ya-1) 


V A 


( ^OX ) inj 0 ' *ox )Ah ox 


+ (1 ^ } inj (1 " *fu J Ah fu + fAp v W) x-0 (^Rtr)^ 


" m ( 5M } 


+ ( *fu ) lnj (1 ' * 


[ (Ap v) x=0 + 

-]] 


; 1nj (1 ’ **>x J 


-L [ 

Av l 


(Ac v) x . 0 ♦ <K ox ) 1nj (1 - *J 


* (S fu> inJ (’ • ♦fu’ j 


» %[' + r^ m t (KR • MR *i ) J 


If the gaseous Injection velocity Is equal to zero (v x , 0 * 0 
state mixture ratio and density at x * Xg are determined by 


fC-46) 


(C -47) 


(C-48) 




(C-53) 

, the steady- 
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(C-51) 


m - m 

*° Vox W 

G \ 1 **„ 

+ Ah fu " m x 0 O 0 ] = m x 0 * 1 


(C-S2) 


These equations were de/eloped by taking the limit as x •* x Q from 
a downstream distance. 


The boundary conditions for the oscillatory differential equations are 
0 x » 0 

P' = Ap (C- 53) 


Fr in these boundary conditions and the oscillatory differential equations 
the oscillatory conditions at the start plane for vaporization (Xg) can be 
determined and are: 






sin 





*0 


..•here 



(C-S4) 


(C-SS) 
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(C-56) 


(C-57) 


(C-58) 


(C-59) 


If the gaseous injection velocity is equal to zero ( v x _q = 0), the oscil- 
latory mixture ratio at x^ is determined by 


MR’ 

*0 




fu illj 0 ♦ m ) - 12. (u.. 7 


A p c . 
x 0 9 


c, ' fu v fu j 

v 


+ 1) 


/m 


A p c. 
x 0 * 


, w; 


vap 


ox 




vap 


ox 


va Pfu 

s — 

vapfu 


- (v; ) AT 


(MR + 1) 
x 0 

(HIT + 2) 

x 0 
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This equation was developed by taking the limit of the mixture ratio 
equation as x -*■ Xq from a downstream distance. 


The nozzle admittance based on upstream conditions is 
\ ■ 


(C-61) 
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and the nozzle admittance based on downstream conditions is calculated 
based on the following analysis. 


The gas flowrate of the nozzle inlet plane is 


. p A t 9 
m = — -x — = Ap v 
c 


where the characteristic velocity is 


c* = 


7 gv RT 


2 -|(y+1)/2(y-1) 


[4 r] 


For short nozzles, the oscillatory mass flowrate can be written as 


m 1 _ p 1 + v 1 _ p' l T 1 / 3c*~ \ MR^ 

m p v p 2 T y » WR" J c* 


Assuming 


£l = 1 El II a I Iii\ El 
p Y P * T \ Y ) p 


the nozzle admittance for a short nozzle can be written as 



(C-62) 


(C-63) 


(C-64) 


(C-65) 


(C-66) 
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Assuming 




MR ■ constant 


the nozzle admittance becomes 




MR' p [ 2r \ 1 

- 1 Pv 


taH IT 


) 


J 



= constant 


(C-67) 


(C-68) 


where Ay can be calculated using the admittance program 

* constant 

developed by Bell (Ref. 32). 
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APPENDIX j 


REPORTS RESULTING FROM CONTRACT NAS9-14315 

As a result of work performed under Contract NAS9-14315, "Orbital Maneuvering 
Engine, Feed System-Coupled Stability Investigation," the following reports have 
been produced. They are listed in chronological order. The subject matter of 
each report is summarized. 

1. Linow, F. R.: "Combustion Characterization Test Plan," ASR 74-333, Rocketdyne 

Division, Rockwell International, Canoga Park, CA, October 1974. Describes 
in detail the test objectives, test hardware, test procedures, and data re- 
quirements necessary to obtain empirical characterization of stability- 
related propellant injection parameters and sensitive operating conditions 
for OME-type combustor hardware. 

2. Kahn, D. R.: "Engineering Model Characterization Evaluation Interim Report," 

ASR 74-372, Rocketdyne Division, Rockwell International, Canoga Park, CA, 
December 1974. A summary of the effort conducted during Phase I of Contract 
NAS9-14315. Includes a detailed assessment of the available techniques for 
modeling the propulsion system's hydrodynamics, combustion dynamics, and 
chamber dynamics, and presents recommendations for the characterization 
technique to be used in each of these areas. A plan for incorporating the 
recommended techniques into an overall engineering model is also described. 

3. Kahn, D. R.: "Engineering Model Verification Plan," ASR 75-108, Rocketdyne 

Division, Rockwell International, Canoga Park, CA, April 1975. Summarizes 
the White Sands Test Facility hot-fire test operating data to be employed in 
the verification of the OME Feed System-Coupled Stability Model. The report 
contains details of Rocketdyne' s thrust chamber and injector, as well as 
schematics of the fuel and oxidizer feedline configuration used during these 
tests. Aoplication of the experimental data to model verification is de- 
scribed and the final verification plan is presented. 
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4. Kahn, D. R., Schuman, M. D., Hunting, J. K., and K. W. Fertig: "Orbital 

Maneuvering Engine, Feed System-Coupled Stability Investigation, Final Report," 
R-980/, Rocketdyne Division, Rockwell International, Canoga Park, CA, Sep- 
tember 1975. This report summarizes all technical effort conducted during 
each phase of the program, and is contained within the present volume. 

5. Schuman, M. D., Fertig, K. W., Hunting, J. K., and D. R. Kahn: "Orbital 
Maneuvering Engine, Feed System-Coupled Stability Investigation, Computer 
User's Manual," R-9808, Rocketdyne Division, Rockwell International, Canoga 
Park, CA, September 1975. Presents complete documentation of the Feed 
System-Coupled Stability Model. The report presents the mathematical formu- 
lation of the model and a detailed description of the computer program. The 
latter includes the structure of the main program and all subroutines, 
instructions for data input, interpretation of program output data, and de- 
tailed analyses of program operation. Appendixes present program flow 
charts, computer code listings, and sample case input/ output data. 
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